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The  meson  component  of  cosmic  rays  in  the  atmosphere  ia  investi- 
gated quantitatively  with  respect  to  its  various  fundamental  charac- 
teristics. The  analysis  is  based  on  a unidimensional  equation  for  the 
vertical  differential  intensity  of  p-mesons,  studied  originally  by 
M.  Sani3.  The  following  topics  are  discussed  in  detail* 

(1)  TK*  a priori  unknown  range  spectrum  of  p-mesons  at  production, 

G(R' ),  is  derived  from  Sands’  equation  with  the  aid  of  recent  experi- 
mental data  concerning  the  p-meson  intensity.  These  include  the  altitude 
and  latitude  dependence  of  p-meson  intensities  as  well  as  the  momentum 

distribution  of  p-raesons  at  sea  level.  It  is  shown  that,  for  the  residual 
-2  -2 

ranges  100  g cm  ^ R'  < 6,000  g cm  , the  production  spectrum  may  be 
represented  by  the  following  empirical  formula* 

G(R' ) ° rg-2c»^sec-18terad-1*] 

(a*R’)3‘58  L 

where  a is  a constant  characteristic  for  a given  geomagnetic  latitude. 

_2 

The  values  of  a vary  from  6U6  g cm  at  the  geomagnetic  equator  to  3 

g cm  at  60  geomagnetic  latitude.  The  value  of  520  g cm  at  50 
geomagnetic  latitude  is  considered  as  the  mo3t  reliable,  since  at  this 
latitude  th9  experimental  data  are  most  complete.  The  production  spec- 
trum ia  compared  with  that  derived  by  M.  Sands  on  the  basis  of  earlier, 
less  accurate  observational  material. 

(2)  The  effects  of  atmospheric  temperature  and  pressure  on  the 
p-meson  intensity  are  studied  for  locations  near  sea  level.  The  treat- 
ment is  rigorous  in  the  sense  that  it  includes  the  continuous  production 
as  well  as  the  ionization  losses  of  p-mescns  in  the  atmosphere.  With  the 
help  of  the  newly  obtained  production  spectrum  and  the  exact  expression 
for  the  survival  probability  of  p-mesons,  a three-terra  regression  formula 
for  the  relative  changes  of  p-meson  intensity  is  derived  and  discussed 

in  detail.  According  to  this  formula,  the  relative  intensity  changes  are 
to  be  correlated  not  only  with  the  average  production  height  and  the  ground 
pressure  (a  customarily  employed  two-term  correlation)  but  also  with  the 
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average  tropospheric  temperature - This  additional  correlation,  resulting 
from  the  ionization  losses  of  p-mesons  in  the  air,  indicates  a possibility 
of  removing  some  apparent  difficulties  in  the  interpretation  of  experi- 
mental data  (as,  for  example,  the  discrepancies  found  in  the  decay  co- 
efficients determined  from  diurnal  and  seasonal  observations,  respectively). 

(3)  The  differential  energy  spectra  of  n-mesons  at  production  are 
computed  for  various  geomagnetic  latitudes,  whereas  use  is  made  of  the 
latitude  dependence  of  the  production  spectrum  of  p -mesons.  The  obtained 
latitude  dependence  of  the  n-meson  spectra  is  linked  with  the  geomagnetic 
effect  on  the  primary  cosmic  radiation.  Some  crude  conclusions  regarding 
the  multiplicity  problem  of  n-meson  production  are  drawn  in  a preliminary 
manner. 
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AN  ANALYSIS  OF  THE  MESON  COMPONENT  OF 
COSMIC  RAYS  IN  THE  ATMOSPHERE 

I.  INTRODUCTION 

Although  as  yet  the  origin  of  the  primary  cosmic  radiation  escapes 
our  knowledge,  we  believe  today  that  vs  understand,  at  least  qualitatively, 
the  secondary  phenomena  caused  in  the  atmosphere  by  the  primary  particle'® 
Let  us  recapitulate  briefly  the  sequence  of  processes  which  are  respon- 
sible for  all  the  complexity  of  the  observed  local  cosmic  radiation. 

The  primary  cosmic  rays  which  consist  of  protons  and,  to  a smaller 
extent,  a-particlea  and  heavier  nuclei,  and  which  are  sufficiently  enei  - 
getio  to  overcome  the  magnetic  field  of  the  earth,  reach  the  upper  at- 
mosphere from  outer  space  and  begin  to  collide  with  the  oxygen  and  nitro- 
gen nuclei  of  the  air.  These  nuclear  collisions  cause  not  only  an  ex- 
citation and  evaporation  of  the  air  nuclei,  but  also  an  emission  of  high- 
energy  secondary  nucleons  as  well  as  a production  of  unstable  particles, 
tt -mesons,  both  charged  and  neutral.*  Neutral  n-mesena,  having  a very 
short  lifetime  (less  than  10*^  sec)  disintegrate  immediately  into  photons, 

and  thus  launch  the  cascade  of  photon-electron  showers.  Charged  n -mesons 

-8 

having  also  a relatively  short  lifetime  (about  2.7  x 10  seo)  decay,  after 
a short  distance  of  travel,  into  neutrinos  and  ^-mesons.  The  latter  par- 
ticles live  on  the  average  2.1  x lO-^  sec  so  that  they  can  travel  a con- 
siderable distance  before  decaying  into  neutrinos  and  electrons.  Further- 
more, they  have  exceedingly  small  cross  sections  for  nuclear  interactions 
so  that  their  energies  can  be  depleted  essentially  only  by  the  ionisation 
losses.  Thus  in  the  lower  parts  of  the  atmosphere  one  expects  the  local 
cosmic  radiation  to  be  composed  primarily  of  two  major  components t the 
electron -photon  component  (also  called  the  soft  component  because  of  its 
absorbability  in  heavy  materials)  and  the  [i -mason  component.  The  third, 
nucleonlo  component,  is  reduced  to  a very  small  fraction  of  its  original 

value  by  nuclear  interactions  (its  absorption  mean  free  path  for  nuclear 

_2 

interactions  is  about  120  g cm  ). 

* We  shall  disregard  here  the  heavier  mesons  since  their  production  rate 
is  very  small  compared  with  that  of  tr -me sons. 
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In  the  following  investigations  we  shall  be  eonoeraeti  with  the 
analysis  of  the  p-meson  component  of  ooemic  rays  in  the  atmosphere. 

This  component a due  to  its  relative  abundance  in  the  lower  atmosphere, 
is  known  in  its  behavior  with  much  greater  accuracy  than  the  remaining 
components  mentioned  above.  Thare  is  a sufficient  amount  of  experi-  * 

mental  material  accumulated  throught  recent  years  to  provide  a reliable 
basis  for  the  mathematical  derivation  of  a quantity  of  primary  importance 
to  cosmic-ray  physics  but  not  accessible  by  direct  measurements  t the  p- 
meson  spectrum  at  production.  In  addition,  there  are  various  observa- 
tional data  which  enable  one  to  test  the  correctness  of  the  assumptions 
made  in  the  calculations. 

However,  it  is  evident  already  from  the  crude  qualitative  picture 
sketched  above  that  the  accurata  mathematical  equations  describing  the 
p-meson  component  must  necessarily  be  of  complex  fora.  It  is,  there- 
fore, almost  impossible  to  arrive  at  usable  results  without  making 
simplifying  assumptions.  Furthermore,  at  some  stages  of  the  formulation 
of  the  problem  one  ia  compelled  to  make  arbitrary  assumptions  because  of 

e 

the  lack  of  satisfactory  theoretical  or  empirical  information.  Briefly, 

even  in  the  oase  of  the  relatively  well-known  p -me son  component  one  muet  * 

be  aware  of  the  highly  limited  validity  of  the  obtained  results.  We 

shall  discuss  this  problem  in  acre  detail  in  Part  II,  Seo.  B. 

We  shall  base  our  analysis  on  a uni  dimensional  model  of  the  dif- 
ferential vertical  intensity  of  p mesons,  the  mathematical  form  of  which 
was  originally  considered  by  M.  Sands  (SMU9).  In  Part  II  of  this  report 
we  shall  derive  the  a priori  unknown  range  spectrum  of  p -mesons  at  produc- 
tion. These  calculations  will  essentially  represent  repetition  of  the 
work  carried  out  by  Sands.  The  necessity  for  re-calculating  the  produc- 
tion spectrum  arises  mainly  from  the  fact  that  the  discrepancies  between 
the  early  experimental  data  used  by  Sands  and  those  obtained  more  recently 
with  higher  accuracy  by  Converai  (CMSO)  and  others  were  too  serious  to  be 
disregarded.  Part  III  deals  with  the  atmospheric  effects  on  the  p-meson 
intensity  near  sea  level.  Using  the  newly  computed  production  spectrum 
we  shall  derive  a regression  formula  for  the  relative  changes  of  the 
p-rceson  intensity  oaused  by  the  variations  of  the  atmospheric  temperature 
and  pressure.  The  results  obtained  will  throw  a new  light  on  the  existing 
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difficulties  In  the  interpretation  of  experimental  data*  Finally,  in 
Parte  17  and  7 we  shall  t email  upon  the  problese  of  the  latitude  dependence 
of  the  |i-  and  n-meson  spectra  at  production.  Some  crude  results  will  be 
utilised  for  speculations  on  the  multiplicity  of  meson  production. 

H.  RANQE  SPECTRUM  OF  p -MESONS  AT  PRODUCTION 

A.  Notation,  Definitions,  Units. 

In  order  to  avoid  repeated  explanation  of  the  notations  used,  we 
list  here  the  symbols  for  the  physioal  quantities  that  will  enter  into 
our  disoussion. 

(a)  Quantities  pertinent  to  the  atmosphere* 

x,s  - atmospheric  depths  at  the  levels  of  production  and  observation 

of  p-mesons,  respectively!  x and  a are  equivalent  to  the  at- 

_2 

mo spheric  pressures  overhead  measured  in  g em  (or  approxi- 
mately in  millibars  }• 

x - atmospheric  depth  at  sea  levelf  x corresponds  to  the  eea-level 


P(x»> 

T(x») 


i (R,b)- 
v 


pressure  and  its  normal  value  is  1,033  go*", 
air  density  as  a function  of  atmospheric  depth, 
temperature  of  the  atmosphere  as  a function  of  atmospheric 
depth,  measured  in  absolute  unite. 

(b)  Quantities  pertinent  to  the  p -mesons 


- 209  a - rest  mass  of  the  p -meson. 

• 

- 2.10  x 10  seo  - mean  life  of  the  p-csson  at  rest. 

- residual  ranges  of  p-mesons  at  the  levels  of  observation  and 

-2 

production,  respectively,  measured  in  g cm  air  equivalent. 

- velocity  of  the  p -me eon. 

- momentum  and  total  energy  of  the  p -meson,  measured  in  multiples 

2 

of  m^c  and  rn^o  respectively. 

- energy  loss  of  p -me son  by  collision  (ionisation  loss)  measured 

2 -2 
in  multiples  of  vs  per  g cm  . 

(o)  Quantities  pertinent  to  the  p -meson  flux* 

)-  differential  vertical  intensity  of  p-mesons|  iT(R,s)dcodR 

represents  the  number  of  p -me sons  which  arrive  at  the 

depth  s from  the  vertical  direction  within  the  solid-angle 

; element,  das,  per  unit  time  and  unit  area,  with  residual 


f • 
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ranges  between  R and  R ♦ dRj  i (R,s)  i*  measured  in  (sec-g- 
sterad)  . Occasionally  the  "short -hand"  tern  "intensity" 
will  be  used  for  iv(R,s). 

IT(Ro,s)-  integral  vertical  intensity  of  p-mesonaj  I (Rq,s)  is  related 
to  iy(R, s)  by  the  equation » cp 
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iv(R,s)dR, 


(1) 


o 
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is,  therefore,  measured  in  (sec-cra  -sterad) 

(d)  Quantities  pertinent  to  the  n -meson t 

276  m - rest  mass  of  the  ir-meson. 

e 

2,65  x 10~°  sec  - mean  life  of  the  w-meson  at  rest, 

momentum  and  total  energy  of  rr -meson  measured  in  multiples 
2 

of  m c and  mo,  respectively, 

TI  TT 

Unidimentional  Equation  for  p, -Meson  Intensity  and  its  Range 
of  Validity. 


Let  us  attempt  to  express  in  mathematical  terms  the  vertical 
p-meson  intensity  which  is  supposed  to  develop  according  to  the  picture 
sketched  in  the  Introduction.  Consider  first  the  number  of  charged 
n-mesons  produced  by  a vertical  flux  of  the  N-component  (mostly  nucleons 
and  n-mesons)  in  an  infinitesimal  layer  dx  at  x with  total  energies  be- 
tween U and  0 ♦ dU  . At  a given  geomagnetic  latitude  this  number  should 
depend  only  on  the  depth  x and  the  energy  U_  since  the  intensity  of  the 

M 

N-component  does  not  change  noticeably  with  time  and  geomagnetic  longitude 
(We  disregard  here  the  abnormal  "storms"  observed  occasionally  in  the  pri- 
mary radiation.)  Now,  It  is  known  from  the  experiments  of  Tinlot  (TJl*8) 
and  others  that  the  component  producing  penetrating  showers  decreases 
exponentially  with  atmospheric  depth  (this  is  true  at  least  for  the  depths 
between  300  and  1,033  S on  ).  Therefore,  if  we  identify  this  component 
with  that  responsible  for  the  rr-meson  production,  we  may  write  for  the 
number  of  tr-mesons  considered  above* 

P(Ufl)e‘*/LdxdUt  , (2) 


where  L is  the  mean  absorption  path  of  the  N-component.  The  function 
P(U  ) may  be  properly  interpreted  as  the  differential  energy  spectrum 

Tf 

of  u -mesons  at  production. 
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As  we  shall  show  below,  ona  may  neglect  the  finite  thickness  of  air 
traversed  by  a w-meson  before  its  decay  into  a p-neson  provided  its  en- 
ergy is  not  exceedingly  high.  By  doing  so  one  can  then  interpret  P(Jff) 
as  the  differential  energy  spectrum  of  tr -mesons  at  decay  wherefrom  one 


oau  obtain  di 


rectly  the  differential  energy  speotrum  of  p-aeaons  at  produc- 


tion M(U).  According  to  Asooli  (AGi>0)  this  latter  quantity  is  related 
to  P(U  ) by: 

If 


where 
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are  the  upper  and  lover  limits  of  the  energies  of  n -masons  that  oan 
give  rise  to  a p-meson  of  energy  U (see  Fig.  16). 

In  what  follows,  it  will  be  more  convenient  to  speak  in  terms  of 
the  differential  range  speotrum  of  p -mesons  at  production,  G(R' ),  rather 
than  in  terms  of  M(U).  These  two  quantities  are  related  by  the  following 
equation: 

G(R' ) - M(U)g,  - M(U)k(U)  . (5) 

One  can  convert  one  opectru*  into  the  other  by  making  use  of  the 
theoretical  curves  for  k(U)  and  the  energy-range  relation  (see  e.g., 

RB$2,  Chapter  II). 

Summarizing,  we  conclude  that  the  number  of  p -mesons  produced  in 

2 

dx  at  x with  residual  ranges  between  R'  and  R'  * dR'  per  sec -cm  -sterad 
may  be  approximated  by: 

G(R')e"X//LdxdR'  . (6) 

In  the  above  estimate  we  have  neglected  the  fact  that  some  of 
the  ir-mesons  as  well  as  some  of  the  p -me sons  will  be  produced  also 
in  directions  different  from  the  vertical.  However,  as  we  shall  see 
below,  for  not  too  low  energies  these  directions  will  differ  only  slightly 


from  the  verticals  Accepting  this  unidimenaional  model  one  notes  that 
at  the  depth  of  observation  the  p-msson  arrives  with  a residual  range 

R - R»  - (s-x)  , (7) 


and  thus  its  differential  vertical  intensity  may  be  written  ast 

8 

iv(R»s)  ■ J G(R+s-x)e  x ^ w(x,s,R)dx,  (8) 

o 

where  w(x, s,R)  represents  the  survival  probability  of  a p -meson  produced 
at  x and  arriving  at  s with  a residual  range  R.  The  explicit  expression 
for  w(x,s,R)  will  be  discussed  in  Sec.  II -D. 

Eq.  (8)  in  its  formal  structure  is  identical  with  that  investigated 
by  M.  Sands  (SMl*9).  Since  the  production  spectrum,  G(R*),  is  not  known 
(either  theoretically  or  experimentally)  Eq.  (8)  represents  an  integral 
equation  which  one  must  solve  with  respect  to  G(R' ) by  making  use  of  the 
observational  data  on  i^(R, a).  Before  we  turn  to  this  problem  let  us  first 
estimate  the  residual  ranges  of  p -me sons  for  which  one  may  consider  Eq.  (8) 
as  a reasonable  approximation  to  reality. 

In  deriving  Eq.  (8)  we  have  made  the  following  assumptions* 

(a)  the  meson  producing  component  varies  exponentially  throughout 
the  entire  atmosphere j 

(b)  the  air  layer  traversed  by  ir-mesons  before  their  decay  is 
infinitesimally  small | 

(c)  the  production  of  w-mesons  and  p -me sons  is  collimated  along 
the  vertlcalf 

(d)  the  multiple  scattering  of  p-mesons  in  the  air  is  negligible. 

Among  all  the  above  assumptions  the  assumption  (a)  is  probably 

the  most  arbitrary  one.  Although  it  is  empirically  justified  for  high- 
energy  events  at  atmospheric  depths  below  the  300-rab  level,  it  is  quite 
uncertain  that  the  same  exponential  law  should  hold  for  lower  energies 
and  for  depths  abov6  300  mb.  We  are  forced  to  make  this  assumption  here 
because  of  the  lack  of  any  better  information.  The  errors  introduced 
thereby  cannot  be  estimated  rigorously.  One  might  only  remark  that  the 
error  should  not  be  too  serious  if  the  significant  deviations  occur  only 
above  100  mb  (this  is  so  because  the  contributions  to  the  p -meson  inten- 
sity originate  mostly  from  the  layers  below  that  level). 


f 
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Concerning  the  assumption  (b)  one  readily  verifies  that  it  repre- 
sents a good  approximation  if  the  momentum  of  e-mesons  is  not  exceedingly 
high.  The  mean  free  path  of  a tr -meson  before  decay  is  given  by* 


Ld  " c^ir?(x)^10’3pirx  (g  cm”2)  (9) 

2 

which  is  of  the  order  of  a few  grams-  per  cm  for  tr-mesons  produced  in 

the  proximity  of  100 -mb  level  with  momentum  smaller  than  100  m c.  For 

tr 

momenta  larger  than  this  value  becomes  comparable  with  the  geometric 

mean  free  path  of  n-raesons  for  nuclear  interaction  and  the  assumption 

that  almost  all  n-mesons  decay  into  p. -mesons  before  interacting  with  air 

nuclei  is  no  longer  valid.  We  shall  take  the  value  of  p *»  100,  which 

-2  n 

corresponds  to  the  n -meson  range  of  about  6,000  g on  , as  the  upper 
limits  for  the  validity  of  Bq.  (8). 

In  contrast  to  the  limitations  on  Bq.  (8)  arising  from  the  assump- 
tion (b)  the  assumptions  (c)  and  (d)  impose  a lower  limit  upon  the  range 
of  validity  of  Bq.  (8).  The  assumption  on  the  collimation  of  tr-mesons 
produced  in  nuclear  interactions  is  essentially  based  on  the  relativistic 
contraction  of  their  cone  of  emission  in  the  laboratory  system.  The  numeri- 
cal value  of  the  maximum  angle  of  emission  with  respect  to  the  line  of 
flight  of  the  primary  particle  depends  on  the  angular  distributions  of 
u -mesons  in  the  oenter-of-mass  system.  Since  the  present  knowledge  on 
this  problem  is  only  of  a speculative  character,  it  is  impossible  to 
make  an  accurate  estimate  of  the  errors  involved.  However,  the  experi- 
mental observations  in  photographic  emulsions  at  high  altitudes  (BRHU9) 
clearly  indicate  that  most  of  the  tr-mesons  with  minimum  ionization  (thin 
secondary  tracks)  are  emitted  in  the  laboratory  system  at  angles  smaller 
than  15°  with  respect  to  the  direction  of  the  primary  particle.  This 
empirical  fact  may  be  considered  as  sufficient  to  support  the  assumption 
on  the  collimation  of  tr-mesons  provided  one  considers  tr-raesons  with  momenta 
larger  than  about  2m^c  (threshold  value  for  minimum  ionization). 

The  maximum  an gle  of  emission  of  p.-mesons  in  the  decay  process  of 
tr-mesons  can  be  readily  calculated  and  is  given  by* 
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Hence  u is  smaller  chan  8°  for  p > 2. 
max  ' 

The  effect  of  the  multiple  scattering  of  p. -mesons  may  be  best 
estimated  from  the  magnitude  of  the  mean-square  angle  which  is  given  by* 
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(E  ■ 21  Mevj  X * 38  g cm”  , see  RB52,  Chapter  II).  The  quantity 
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appearing  in  the  denominator  of  the  integrand  in  Eq.  (11)  is  practically 

-2-1  2v 

constant  for  all  p ,>  2 (its  numerical  value  for  air  is  2.1  x 10  g cm  ) 
so  that  one  has  the  following  simple  formula  for  the  mean-square  angle: 


<3,  * °-o5°  ( m~  m;  ) • 


Eq.  (12)  implies,  for  example,  that  a p -meson  produced  at  100-mb  level  and 
arriving  at  sea  level  with  p * 3 will  have  a rms  angle  of  7.0°.  For  momenta 
smaller  than  this  value  the  deviations  from  the  vertical  become  more 
critical. 

Summing  up  we  see  that,  if  we  exclude  p -mesons  with  momenta  smaller 
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that  3 m^c  (corresponding  range  100  g cm  ) we  may  expect  that  the  local 

p-meson  intensity,  produced  by  a vertical  flux  of  the  N-component,  will 

be  contained  within  a vertical  cone  of  30°  zenith  angle.  The  uni- 

dimensional  treatment  of  the  development  of  the  p-meson  component  can 

be  considered  valid  within  these  limits. 

Referring  to  the  latter  limitation  and  to  the  limitation  arising 

from  the  assumption  (b),  we  conclude  that  the  range  interval  in  which 

Eq.  (8)  represents  a reasonable  approximation  to  reality  extends  from 
-2  -2 

R ■=  100  g cm  up  to  R ■ 6,000  g cm  air  equivalent. 
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C.  'Experimental  Information  Concerning  the  ^i-Keeon  Component . 

As  we  have  mentioned  in  the  foregoing  section,  Bq.  (8)  can  be  solved 
with  respect  to  the  range  spectrum  at  production  only  if  one  has  at  his 
disposal  a sufficient  amount  of  experimental  data  concerning  the  p -meson 
intensity,  iv(R,a).  We  have  chosen  the  following  two  experiments  as  the 
most  reliable  and  suitable  for  the  purposes  of  our  analysis. 

(a)  measurements  of  the  differential  di stribution-in -momentum 
of  cosmic-ray  mesons  at  sea  level  by  Caro,  Parry  and  Rathgeber  (CDR50)j 

(b)  measurements  of  the  altitude  dependence  of  cosmic-ray 

_2 

mesons  with  rssidual  ranges  between  100  and  117  g cm  air  equivalent 
by  M.  Conversi  (CM$0). 

In  the  first  experiment  the  penetrating  cosmic-ray  particles  coming 
nearly  from  the  vertical  direction  (the  opening  angle  was  about  16°)  were 
deflected  in  the  air  gap  of  an  electromagnet  and  the  deflection  was  re- 
corded by  means  of  an  array  of  G.M.  counters.  The  measurements  were  car- 

2 

rled  out  at  sea  level  and  covered  the  momentum  range  between  2.U  x 10 
and  5 x 10^  Mev/c.  The  differential  momentum  spectrum  thus  obtained  is 
reproduced  in  Fig.  1.  Since  the  particles  recorded  represent  almost  ex- 
clusively the  vertical  p. -me  sons,  one  can  utilize  the  above  data  for  re- 
constructing the  differential  vertical  intensity  of  p-mesons  at  sea  level, 
iv(R,xQ),  as  a function  of  the  residual  range  R.  Applying  the  theoretical 
expression  for  the  momentum  loss,  k/p,  for  |i-mesons  in  air  we  arrived  at 
a curve  shown  in  Fig.  2.  The  solid  line  corresponds  to  the  measurements 
by  Caro  et  al,  while  the  dashed  line  represents  Rossi's  curve  (RBl|8)  drawn 

for  comparison.  Note  significant  deviations  from  Rossi's  curve  at  residual 

-2 

ranges  above  3,000  g cm  , 

In  addition  to  the  above  sea-level  data,  we  shall  find  the  measure- 
ments by  Conversi  as  indispensable  for  studies  concerning  the  behavior 
of  the  production  spectrum  at  small  residual  ranges.  Combining  the 
techniques  of  delayed  coincidences  and  anticoincidences,  Conversi  was 
able  to  determine  the  absolute  number  (per  sec-g-sterad)  of  p.-mesons 
stopped  in  10  cm  of  graphite  after  traversing  15.2  cm  of  lead  for  various 
atmospheric  depths  at  50°  N geomagnetic  latitude.  The  counter  telescope 
had  a relatively  good  geometry  (the  maximum  zenith  angle  permitted  for  a 
meson  entering  the  telescope  and  stopping  in  the  absorber  was  about  35°) 
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figure  1 

Differential  momentum  distribution  of  u-mesons 
at  sea  level  measured  by  Caro  at  al.  (CDiftO). 
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and  the  statistical  accuracy  of  the  data  exceeded  that  of  earlier  measure- 
ments of  this  kind  (RBU?>  SMli9)  considerably.  The  delay  curves  yielded 
a mean  life  which  was  consistent  with  ohat  of  y -mesons  at  all  depths  con- 
sidered, Thus  one  may  interpret  Conversi's  data  as  a direct  measure  for 

the  dependence  of  the  vertical  differential  intensity  on  the  atmospheric 
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depth  for  a y -meson  with  residual  ranges  between  100  and  117  g cm  air 
equivalent „ Conversi's  results  with  corresponding  errors  are  shown  in 
Fig,  3f  the  solid  line  represents  the  differential  vertical  intensity  as 
computed  on  the  basis  of  the  production  spectrum  derived  in  Sec,  II-E. 

The  dashed  line  represents  the  normalized  intensity-depth  curve  deduced 
by  Sands  (SMU9)  from  his  earlier  measurements.  The  large  difference  in 
the  behavior  of  Sands'  and  Conversi's  curves,  respectively,  is  primarily 
due  to  the  fact  that  Sands'  data  cannot  be  interpreted  as  directly  repre- 
sentative for  the  vertical  differential  intensity  of  y -mesons.  They  rather 
represent  an  intensity  integrated  not  only  over  all  ranges  of  mesons  be- 
tween  $ and  83  g cm  but  also  over  all  directions  of  incidence  (the  opening 
angle  of  Sands'  arrangement  was  dose  to  180°), 

Concluding,  we  would  like  to  call  the  reader's  attention  to  the  good 

agreement  existing  betwoen  the  measurements  of  Caro  et  al.  and  those  of 

-2 

Conversi  for  y-raesons  at  sea  level  with  residual  ranges  of  100  g cm  , 

(Note  that  both  experiments  were  performed  at  geomagnetic  latitudes  near 
5 0 ° although  on  opposite  sides  of  the  geomagnetic  equator.) 

D.  Survival  Probability  of  y -Mesons . 

Another  quantity  which  enters  into  our  basic  equation  (8)  is  the 
survival  probability  of  y-mesons,  w(x, s,R).  Due  to  its  implicit  de- 
pendence on  the  vertical  distribution  of  the  atmospheric  temperature, 
this  quantity  will  be  treated  separately  in  some  detail. 

For  a y-meson,  produced  at  the  elevation  a,  which  travels  verti- 
cally toward  the  earth  with  the  velocity  cp  and  is  observed  at  the  eleva- 
tion zq  the  survival  probability  Is  given  byi 
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Figure  3 

Vertical  differential  intensity  of 
p-mesons  with  residual  ranges  be- 
tween 100  and  11?  g cm"-  air  equiv- 
alent as  a function  of  atmospheric 
depth  at  50°  N geomagnetic  latitude. 
The  full  circles,  with  corresponding 
errors,  represent  the  measurements 
of  Con\ersi  (CM£0)j  the  open  circle 
represents  the  measurements  of  Xrau- 
shaar  Q(WU9)  at  Fcho  Lake.  The 
dashed  line  represents  Sands'  mea- 
surements on  p -mesons  with  residual 
ranges  between  5 and  83  g cm"2.  The 
solid  curve  was  computed  by  means  of 
Eqs.  (11-32)  and  (11-33). 
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the  factor  Tl-p2  in  Eq.  (13)  accounts  for  the  relativistic  time  dilata- 
tion. It  is  conventient  to  express  the  distance  element,  dz',  in  terms 
of  the  increment,  dx' , of  the  atmospheric  depth.  Since  dx’  ■ ^(x^dz’ 
Eq.  (13)  may  be  written  as* 


w(x,s,R)  “ exp  J J 


; i 


Here  p'  stands  for  the  momentum  of  the  p, -meson  at  the  depth  x'  expressed 
as  a function  of  the  corresponding  residual  range  (R  ♦ s - x' ). 

We  shall  find  it  useful  for  the  consideration  in  Part  III  to  ap- 
proximate the  momentum-range  relation  by  the  following  analytical  formula* 


1 - B - U 

p FTT?  ^ * 


-2  -2  -2 
where  B ■ 53,5  g cm  j b ■ 56  g cm  if  R is  measured  in  g cm  air 

equivalent;  J6  ■ constant  - 2.07  x 10"3.  With  the  numerical  valueB  of  the 

constants  quoted  above,  Eq.  (15)  is  applicable  for  all  pi -mesons  with  resi- 
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dual  ranges  between  30  and  6,000  g cm  air  equivalent.  In  this  region 
it  reproduces  the  theoretical  curve  (HB52,  pageB  1*0— 1*1)  within  an  accuracy 
of  one  percent.  Referring  to  the  discussion  of  the  validity  of  Eq.  (8)  in 
Sec.  II-B  we  may  consider  Eq.  (15)  as  sufficient  for  the  purposes  of  our 
analysis. 

The  evaluation  of  Eq.  (li*)  also  requires  the  knowledge  of  the 
vertical  behavior  of  the  air  density,  <p(x*)  or  of  the  atmospheric  tem- 
perature, T(x').  These  two  quantities  are  related  by  the  following 
equation* 

*'  _ ^T(l') 

-pF)  T <l6) 

where  ^ is  the  universal  gas  constant,  M is  the  mean  iuoio  c uXs.  r weight 
of  air,  and  g is  the  acceleration  of  gravity  ( $,/Mg  - 2.87  x 103  cm/°C). 
It  is  more  appropriate  to  discuss  the  survival  probability  in  terms  of 
T(x' ) rather  than  in  terms  of  ^)(x' ) because  T(x' ) i3  a quantity  that 
can  be  measured  directly  by  means  of  radiosondes. 

By  combining  Eqs.  (li*),  (15)  and  (16)  one  can  write  for  the 
survival  probability* 
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* > 
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,U,e,R)  - «p  [-  J Il|i  (JL^  -\t)  fcj], 


where  we  have  put  for  short t 


x » R ♦ s ♦ b. 
0 


Since t 


11  *,  * v**  J 


the  following  representation  of  w(x,s,R)  is  possible! 

w(x,s,R)  - exp |odR( s, R)H(x,e)  + OCK(s,R)K(x,s,R)Jjf 


where 


°V'-R>  ■%  (\  -<«)■-  orp(.^)  ’ 


(17) 


(18) 


(19) 


OiE(  s,R)  - - S g , 

H MgcTx/ 


(20) 


and 


H(x,s)  ( ^)dx«  , 

MgJ 


(21) 


K(x,s,R)  - xm 


( m~- 


(22) 


The  representation  of  x(x,s,R)  given  by  Eq.  (18)  has  some  advantages 
over  those  given  by  Eqs.  (lit)  or  (17).  The  two  terms  in  the  exponent 
of  Eq.  (18)  have  a direct  physical  significance.  The  function  H(x,s) 
defined  by  Eq.  (21)  represents  simply  the  distance  from  the  pressure 
level  x (tho  production  level)  to  the  pressure  level  s (the  observation 
level) . 
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The  first  term  in  Bq.  (18),  exp(CX^H),  is  the  main  term;  it  re- 
presents the  survival  probability  for  a ^ -meson  produced  at  x and  travel- 
ing the  distance  H with  a fixed  momentum  p(s  ♦ R).  For  most  of  the  mesons 
recorded  in  the  lower  atmosphere,  x is  considerably  smaller  than  s so  that 
p(s  ♦ R)  will  not  differ  significantly  from  the  actual  momentum  at  produc- 
tion p(s  + R - x).  This  implies  that  approximately 


” “P  ['  • <23> 

One  recognizes  in  the  right-hand  side  expression  a formula  that  is  often 
quoted  in  the  literature.  It  represents  the  survival  probability  of 
tx-mesons  if  one  neglects  their  ionization  losses  in  the  air. 

The  second  term  in  Eq.  (18),  exp (06 ^K),  may  be  regarded  as  a cor- 
rective terra  accounting  for  the  ionization  losses  of  (x-mesons  in  the  air. 
The  function  K(x,s,R)  defined  by  Bq.  (22)  will,  in  general,  depend  not 
only  on  the  temperature  distribution  but  also  on  the  residual  range  of 

the  p. -meson  -under  consideration.  How«n •.  -,r.  for  sufficiently  large  R 

m,o 

(say,  R >1,000  g cm  ) K(x,s,R)  may  be  written  roughly  as 


K(x,s) 


T(x')dx» 


x 


x'dz  , 


(2U) 


i.e.,  in  case  of  fast  mesons,  the  function  K is  closely  related  to  the 
temperature  averaged  over  the  region  between  the  levels  of  production 
and  observation  or,  what  is  equivalent,  to  the  amount  of  air  mass 
overhead. 

For  the  numerical  evaluation  of  Eq.  (18)  one  needs  the  meteorologi- 
cal data  on  T(x’ ) which,  in  general,  will  be  different  for  different 
seasons,  different  geographic  locations,  etc.  However,  a closer  study  of 
meteorological  data  shows  that  for  the  description  of  the  annual  means 
of  atmospheric  temperature  one  needs  essentially  only  two  parameters t 
the  atmospheric  depth  (i.e.,  the  pressure  overhead)  and  the  geographic 
latitude.  The  longitudinal  variations  of  the  annual  means  at  a given 
latitude  circle  turn  out  to  be  insignificant  except  for  the  localities 
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with  extreme  weather  conditions  m,  e-g.,  localities  with  a pronounced 


continental  or  maritime  climate.  Eren  in  those  exceptional  oases  the 
deviations  become  significant  only  in  the  lower  parts  of  the  troposphere. 
Hence,  it  is  plausible  to  oonsider  the  annual  means  of  T(x' ) averaged 
along  a given  latitude  circle  as  basic  quantities.  The  local  corrections 
may  be  then  oarried  out  with  the  aid  of  a regression  formula.  An  explicit 
expression  of  such  a formula  will  be  discussed  at  length  in  Part  III  of 
this  report.  Here  we  shall  be  merely  concerned  with  the  annual  means  of 
T(x')«  Fig*  1*  shows  the  behavior  of  these  means  plotted  as  functions  of 
the  atmospheric  depth  for  various  geographic  latitudes.  The  curves  were 
constructed  from  a graph  given  by  Haurwitz  and  Austin  in  their  textbook 
on  climatology  (HBlUi,  page  37).  Note  the  inverse  behavior  of  the  tem- 
perature in  the  stratosphere  and  the  troposphere  when  one  passes  from  the 
equatorial  region  to  the  middle  latitudes.  The  temperature  distribution 
in  the  upper  stratosphere  (above  10  mb)  is  not  as  yet  definitely  known; 
however,  this  uncertainty  does  not  affect  our  considerations  noticeably. 

By  means  of  Eq.  (18)  and  Fig.  U one  can  calculate  numerically  the 
average  survival  probability  for  any  given  latitude.  Since  the  experi- 
ments of  Conversi  and  Caro  et  al.  were  oarried  out  in  the  proximity  of 
U0°  geographic  latitude,  we  based  our  calculations  on  the  temperature 
distribution  of  this  region.  Referring  to  the  experiment  of  Conversi 
we  show  in  Fig.  5 the  survival  probability  for  n. -mesons  produced  at  the 
depth  x and  observed  at  the  depth  s with  residual  ranges  of  100  g cm 
Q*(x, s,100)  is  plotted  versus  the  amount  of  air  traversed,  s-x}.  Refer- 
ring to  the  experiment  of  Caro  at  al.  we  show  in  Fig.  6 the  negative 
natural  logarithm  of  survival  probability  for  n -mesons  produced  at  x and 
reaching  sea  level  with  various  residual  ranges  between  100  and  2,000  g cm  . 

Regarding  the  n -mesons  which  arrive  at  sea  level  with  residual  ranges  lar- 
_2 

ger  than  2,000  g cm  , it  is  useful  to  approximate  their  survival  probab- 
ility by  the  following  simple  formula t 


w(x,xq,R) 


v(R) 


(25) 


o 

where  the  quantities  v(R)  and  ^(R)  are  numerically  derived  functions 
of  the  residual  ranges  at  sea  level  and  are  shown  in  Fig.  7.  The  right- 
hand  side  expression  in  Eq.  (25)  reproduces  the  survival  probability 
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very  accurately  for  all  production  levels  with  x ^ 500  g era  j at  pro- 
duction levels  closer  to  the  ground  (contributing  very  little  to  the  value 
of  Eq.  (8)  ) the  formula  (25)  introduces  errors  which  amount  only  to  a few 
percent  of  the  actual  values  of  w(x,x  ,r'„  It  is  interesting  to  compare 
Eq.  (25)  with  the  survival  probability  of  fast  mesons  in  the  case  of  the 
isothermal  atmosphere.  The  latter  probability  may  be  obtained  directly 
from  Eqs.  (21)  and  (23)  by  setting  T ■ constantj  them 


x<fcT/cl:Mgp 

wisothermal  ^ 

o 


(25*) 


The  resemblance  between  Eqs.  (25')  and  (25)  is  emphasized  by  the  fact 
that  the  function  v(R)  does  not  depart  markedly  from  unity,  and  the 
function  (R)  behaves  similarly  to  the  expression  $T/c'CMgp(xo+R). 

E.  Solution  of  the  Intensity  Equation  with  Respect  to  the  Produc- 
tion Spectrum. 

We  are  now  prepared  to  turn  to  our  main  problem t the  solution  of 
Eq.  (8)  with  respect  to  the  range  spectrum  of  -mesons  at  production. 
Following  the  procedure  of  Sands  (SHh9)  we  find  it  convenient  to  begin 
with  the  determination  of  G(R! ) at  large  values  of  residual  ranges. 

(a)  Production  spectrum  at  large  residual  ranges. 

If  one  considers  the  \l -meson  intensity  only  for  large  ranges  at  sea 

-2% 

level  (say,  R)  2,000  g cm  ) one  can  introduce  two  simplifications  in 
Eq.  (8)j  first,  one  can  express  the  survival  probability  in  terras  of 
Eq.  (25),  second,  one  can  assume  the  production  spectrum,  G(R’ ),  to  be  a 
slowly  varying  function  within  the  interval  of  integration.  The  latter 
assumption  is  justified  by  the  fact  that,  for  large  R,  the  residual  range 
at  production,  R1  ■ R+xq-x,  does  not  change  strongly  within  the  interval 
0 ^x^  xq.  The  slow  variation  of  G(R’ ) suggests  the  application  of  the 
mean  value  theorem  of  integral  calculus  to  the  right-hand  side  of  Eq.  (8). 
According  to  the  discussion  in  the  Appendix  one  has  then 

xq 

iy(R,xo)  - G(R+xQ-xm)  J e_X//L  w(x,xo,R)dx.  (26) 

o 

The  substitution  of  Eq.  (25)  yields  for  Eq.  (26) 


iy(R,xo)  - G(R+xo-xm)Lw(L,xo,R)T(l+  /\  /x/L)  > 


(27) 


f j .**■  ■ 


» 


'V* 


where  y 

Ttay)  - y"  tz“1e"tdt  T(a)  “ y*”1®"7  (28) 

is  the  incomplete  gamma  function.  Due  to  the  fact  that  xq/L  is  con- 
siderably larger  than  one,  one  can  approximate  T"*(l+^  jxq/L)  by  the 
complete  gamma  function,  V (1+  ).  The  value  of  the  "mean*  depth 

x can  be  estimated  with  the  aid  of  Eq.  (A-2)$  one  Lnds  thus* 


T1;  2*^j  x /l) 

x»-  t&T,  x°A)  (1,?>)L  • (Z9) 

Hence,  referring  to  Eqs.  (27)  and  (29),  one  has  the  following  ap- 
proximate solution  for  the  production  spectrum  at  large  ranges  t 


Q(R+x  -x  ) 
o m 


VR*X0) 

"iw(t,xo,li)T'0>T)  * 


(30) 


Taking  for  the  absolution  mean  free  path,  t,  ths  numerical  value 
-2 

of  121  g cm  and  using  Figures. 2 and  7,  we  have  computed  G(R+x  -x  ) 

-2  —2  ® ® 
for  2,000  g cm  < R < 6,000  g cm  . The  resulting  curve  is  shown 

in  Fig.  8.  At  first  sight  this  curve  seems  to  be  expressible  in  terras 

of  a simple  power  law.  However,  a closer  inspection  shows  that  the 

following  representation  is  more  accurate* 


Q(R»)  - (3D 

(a+R)' 

where  C,  a and  y are  positive  constants.  Unfortunately,  the  range 
interval  for  which  Eq.  (30)  is  valid  is  too  small  to  permit  an  unam- 
biguous choice  of  a unique  set  of  the  three  constants  C,  a,  and  ^ . 

One  can  readily  verify  that  it  is  possible  to  reproduce  the  curve 

_2 

in  Fig.  8 by  assuming  any  arbitrary  value  between  200  and  700  g cm 
for  the  constant  a,  if  the  remaining  constants,  C and  , are  properly 
adjusted.  Table  I gives  four  such  s?ts  of  the  numerical  values  of  C, 
a,  and  ^ , which  reproduce  the  spectrum  curve  in  Fig.  8 equally  well. 
Of  course,  any  set  of  constants  obtained  by  an  interpolation  of  the 
values  quoted  in  Table  I is  also  compatible  with  our  solution  Eq.  (30). 
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Eange  spec  trim  of  ^x-naaona  at  production  for 
large  residual  ranges. 


TABLE  I 


Four  sets  of  constants,  a,  ^ , and  C,  which  are  equivalent  for  the 

empirical  representation  of  the  production  spectrum  of  (i-mesons  at 

_o 

residual  ranges  larger  than  3*000  g cm  air  equivalent: 

G(R' ) - C(a+R')“*  „ 


Set 

No. 

a 

t ~2n 

(g  cm  ) 

r 

G 

(g*~2cra2“2/  sec^sterad"1) 

1. 

300 

3.1:6 

2.25  x 101* 

2. 

Uoo 

3.51 

3.30  x 101* 

3. 

5'00 

3.57 

6.61  x 10^ 

U. 

600 

3.63 

11.8U  x 10^ 

'*.*S**W*M>  • V 
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In  order  to  eliminate  this  ambiguity  one  must  turn  to  the  region  where 

the  residual  range  at  production  assumes  values  comparable  with  those 

_2 

permitted  for  the  constant  a,  i.e.,  smaller  than  1,000  g cm  , Evi- 
dently, the  p -meson  intensities  measured  by  Conversi  correspond  to  such 
a region.  Thus,  by  making  use  of  Conversi' s data  one  can  test  if  the  pro 
duction  spectrum  at  small  ranges  is  expressible  by  a formula  suggested 
from  the  large-range  data  and,  if  so,  to  determine  a unique  set  of  values 
for  the  characteristic  constants. 

(b)  Production  spectrum  at  small  residual  ranges. 

According  to  the  discussion  in  the  foiegoing  sections  we  may  trans- 
late the  data  of  Conversi  into  the  following  equation t 

a 

iy(l00,s)  - y*  0(100+s-x)  e _3^L  w(x,s,100)dx  ? (32) 

o 

where  the  survival  probability  is  given  by  Fig.  *>. 

Assuming  that  the  production  spectrum,  0(100*3 -x),  is  expressible 

by  Eq.  (31)  we  have  evaluated  numerically  the  right-hand  side  of  Lq.  (32) 

_2 

for  various  depths  s ranging  between  200  and  1,000  g cm  . The  calcula- 
tions were  repeated  four  tiraos,  each  time  with  a different  set  of  con- 
stants a,  C,  and  ^ quoted  in  Table  I,  The  results  thus  obtained  were 
then  compared  with  the  experimental,  points  given  in  Fig,  3.  Among  the 
computed  curves,  the  one  corresponding  to  the  third  set  of  constants 

r I 

(a  ■ !>00  g cm  J ■ 3.57,  C * 6.6i  x 104)  was  found  to  be  in  closest 
agreement  with  the  experimental  data.  When  a set  of  constants  obtained 
by  an  interpolation  between  the  third  and  fourth  sets  of  Table  I was 
taken,  the  agreement  became  practically  complete  except  for  one  point 
at  the  depth  of  231  g era  (see  the  solid  line  in  Fig*  3).-.  Conversi' s 
value  of  the  p -meson  intensity  at  this  depth  falls  above  the  calculated 
curve.  We  believe,  however,  that  this  measurement  should  be  corrected 
for  the  p -mesons  locally  produced  in  the  lead  absorber  by  the  nucleons 
which  are  still  abundant  at  this  high  altitude  (for  remarks  on  this 
problem  see  also  CM50).  Corrections  of  this  kind  may  be  expected  to  be 
of  the  order  of  20  percent.  Therefore,  we  do  not  consider  this  single 
discrepancy  as  crucial. 


-2.7* 


Summarizing,  we  conclude  that  the  following  expression  for  the 
p-raeson  spectrum  at  production* 


G(R>) 


7.31  x IQ1* 

(520  ♦ R»)3*58 


2 

era  sec 


-1 


aterad 


-il 


(33) 


is  oonpatible  with  the  experimental  data  for  residual  ranges  between 

-2  -2 
100  and  1,100  g cm  (Converai's  data),  and  3»000  and  6,000  g cm 

(Caro's  data).  In  order  to  test  the  compatibility  of  the  above  formula 

-2 

also  for  the  intermediate  ranges  between  1,000  and  3,000  g ora  it  is 

sufficient  to  calculate  the  values  of  the  differential  intensity  of 

p -mesons  which  arrive  at  sea  level  with  residual  ranges  between  100  and 
_2 

2,000  g ora  and  to  corapare  the  results  with  Fig.  2.  This  has  been  done 

numerically  by  using  the  graphs  in  Fig.  6 for  the  survival  probability. 

The  computed  points  are  indicated  by  open  oircles  in  Fig.  2.  One  sees 

that  the  agreement  is  fully  satisfactory. 

In  Fig.  9 we  show  the  p-raeson  spectrum  at  production  computed  from 

Eq.  (33)  for  small  and  intermediate  ranges*  the  dashed  curve  represents 

the  production  spectrum  as  dorived  by  Sands.  It  is  interesting  to  note 

that  our  smooth  production  spectrum  behaves  in  such  a way  as  to  average 

-2 

out  the  "dip"  obtained  by  Sands  at  R’  ■ 100  g era  . 

We  close  the  discussion  of  this  part  cf  our  analysis  with  a few 

remarks  concerning  the  degree  of  accuracy  that  cne  may  claim  for 

—2 

Eq.  (33).  Since,  for  2,000  < R<  6,000  g cm  , Eq.  (33)  is  equivalent 
to  Eq.  (30)  we  can  estimate  the  relative  error  made  in  our  computation 
of  G at  large  ranges  by  means  of  Eq.  (A-U),  viz* 


1 2 
relative  error  ■ ^(l+^jL 


2 

1 dG 


dR« 


R' 


R+x  -x 
o m 


- ) (1+*) 


a+R+x  -x 
o m 


— 

in_ 


2 
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For  R > 2,000  g cm  this  error  is  smaller  than  1.2  percent.  Hence, 
referring  to  Eq.  (30),  at  large  residual  ranges  one  may  consider  the  pro- 
duction spectrum  to  be  known  almost  as  accurately  as  the  p -meson  intensity 
at  sea  level  itself.  It  is  more  difficult  to  estimate  the  errors  contained 
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Bange  spectrum  ofU-meBon#  at  production  for  «mall  and  inter- 
mediate residual  ranges.  The  daehed  curre  represent*  the  pro- 
duction spectrum  derired  hy  Sand#  (SM49).  The  solid  curve  wa* 
computed  from  Xq.  (11-33) 
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in  the  determination  of  G at  small  ranges.  However,  the  accuracy  with 
which  the  constant  a is  known  can  be  taken  as  a good  measure  for  the 
uncertainties  involved.  Our  repeated  evaluations  of  Eq.  (32)  for  dif- 
ferent values  of  a showed  that  1^(100, s)  is  quite  sensitive  to  the  as sinned 

values  of  this  constant.  For  example,  a 5 percent  uncertainty  in  i (100, s) 

-2  v 
at  s * 300  g cm  produced  only  a 2 percent  uncertainty  in  a.  Since  the 

exponent  ^ is  relatively  insensitive  to  the  values  of  a (see  Table  I)  and 
is  essentially  determined  by  large-range  measurements,  one  can  regard  ex- 
pression (33)  for  the  production  spectrum  at  50°  geomagnetic  latitude  as 
comparatively  well  established. 

IH.  ATMOSPHERIC  EFFECTS  ON  p-MESON  INTENSITY  NEAR  SEA  LEVEL 

A.  Oeneral  Remarks. 

Numerous  experiments*  of  the  last  two  decades  have  shown  that  the 
intensity  of  cosmic  rays  is  influenced  by  the  atmospheric  conditions 
existing  during  the  period  of  observation.  It  has  been  found  from  ex- 
periments on  the  ground  that  the  variations  of  the  cosmic -ray  intensity 
are  closely  related  to  the  changes  of  pressure  and  temperature  of  the 
atmosphere  above  the  observer.  While  the  pressure  effect  simply  indicates 
the  dependence  of  cosmic  radiation  on  the  amount  of  air  mass  traversed, 
the  temperature  effect  has  its  origin  in  the  instability  of  p -mesons 
(BFM38)  whioh  form  the  preponderant  part  of  the  penetrating  component  of 
cosmic  rays  at  sea  level.  The  following  simple  argument  shows  qualitatively 
that  an  increase  of  the  temperature  causes  a decrease  of  ths  p-meson  in- 
tensity at  ground  level.  Each  p-raeson  produced  at  a certain  atmospheric 
depth  - i.e.,  at  a given  pressure  level  - has  to  travel  the  distance  from 
the  production  layer  to  the  ground  in  order  to  be  deteoted.  An  increase 
of  the  temperature  will  increase  this  distance,  and  thus  enhance  the  pro- 
bability of  decay  in  flight.  If  all  mesons  were  produced  at  the  same  at- 
mospheric depth,  and  if  one  could  neglect  their  ionization  losses,  the 
fluctuations  in  the  height,  H,  of  a single  pressure  level  - the  production 
level  - would  suffice  to  account  for  the  temperature  effect  on  the  p -meson 
intensity.  One  could  then  express  the  variations  of  the  integral  p-meson 

* For  a review,  see,  e.g.,  H.  Elliot,  Chapter  VIII  in  "Progress  in  Cosmic 
Ray  Physics"  ( North-Holland  Publishing  Company,  1952). 
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intensity,  I,  with  an  equation  of  the  following  type: 

Sl/I  - <f  H + ^xo 

where  <$*  H and  Oxq  a^  the  variations  of  the  height  of  production  and 
of  the  ground  pressure,  respectively.  The  coefficients  A^  and  are 
often  referred  to  as  the  "decay”  and  "absorption”  coefficients,  respectively. 

Early  investigations  seemed  to  confirm  the  approximate  validity  of 

the  above  regression  formula.  A.  Duperier  (DAUU,  DAl+8)  in  particular, 

found,  through  the  statistical  analysis  of  the  observed  data,  that  the 

ground  pressure  and  the  height  corresponding  to  the  atmospheric  depth 
_2 

of  10U  g cm  were  the  controlling  factors  in  the  variations  of  cosmic- 
ray  intensity  at  sea  level.  However,  the  more  recent  studies  by  the  same 
author  (DAh9),  as  well  as  those  by  D.W.H.  Dolbear  and  H.  Elliot  (DDW?1), 
have  shown  that  the  two-term  equation  (1)  is  inadequate  to  account  fully 
for  the  variations  of  the  cosmic -ray  intensity  at  sea  level.  The  incon- 
sistencies have  not  been  explained  quantitatively  in  a satisfactory  manner. 

In  what  follows,  we  shall  attempt  to  give  a more  rigorous,  quantitative 
treatment  of  the  problem  outlined  above.  In  particular,  by  making  use  of 
the  intensity  equation  (II-8)  we  shall  be  able  to  take  into  account  two 
facts  customarily  neglected  in  rough  estimates  of  atmospheric  effects: 

(l)  the  fact  that  the  n-mesons  are  produced  continuously  throughout  the 
atmosphere,  and  (2)  the  fact  that  the  n -mesons  suffer  ionization  losses 
during  their  propagation  through  air.  As  we  shall  see  in  the  following 
sections,  the  first  fact  will  justify  and  clarify  the  notion  of  an  average 
production  layer  for  the  bulk  of  p-mesons,  and  the  second  fact  will  intro- 
duce a third  term  into  the  regression  formula  discussed  above.  This  addi- 
tional term  will  turn  out  to  be  primarily  controlled  by  the  temperature 
changes  in  the  lower  atmosphere,  where  the  ionization  losses  are  relatively 
more  important  than  in  the  upper  atmosphere. 

All  considerations  in  this  part  apply  only  to  the  vertical  intensities 
(both  differential  and  integral)  of  n-mesons  measured  at  locations  near 
sea  level. 
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B.  Ten 


Effects  on  n -Me a on  Intensity. 


Let  us  discuss  first  the  case  where  the  ground  pressure,  xq,  is 
kept  constant.  Then  the  variation  of  i^(R,XQ)  due  to  the  temperature 
changes,  8 T(x' ),  is  given,  according  to  Bq.  (II-8)  by» 


XQ 

6Tiv.(R,xQ)  - J'  ((fTv/w)  $ (x,xo,R)dx, 


where 


r -oc/L 

<P  (x,xq,R)  - G(R+xo-x)e  w(x,xq,R), 

and,  according  to  Eq.  ( 11-18), 

S^/w  - OdR(xo,R)  $ H(x,xq)  + OCk(xo,R)  <5"k(x,xo,R), 

V***.)  • f°  %-)d"  • 


&K( x,x  ,R)  ■ 


(xQ+R+b)  J' 
x 


otix' ; 

x ♦R+b-x' 

o 


d ’ 


In  order  to  evaluate  Eq.  (2)  we  shall  make  use  of  the  mean  value  theorem 

discussed  in  detail  in  the  Appendix. 

Evidently,  the  function  , defined  by  Eq.  (21 ),  satisfies  the  con- 

_2 

dition  (2)  of  the  Appendix  very  well  for  R >100  g cm  ; furthermore, 
according  to  Eq.  (JI-22),  the  function  5* X satisfies  the  condition  (1). 
Hence,  we  may  apply  the  mean  value  theorem  directly  to  the  second  term 
in  Eq.  (2),  We  cannot  apply  it  directly  to  the  first  term,  because  the 
function  5 H displays  a logarithmic  divergence  at  x ■ 0.  However,  if 
we  consider  separately  the  bounded  functions  8 H/ln(x  /x)  and  & ln(x  /x), 
rather  than  the  functions  and  , we  readily  verify  that  these  two 


t 
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functions  do  satisfy  conditions  (1)  and  (2),  respectively.  Consequently, 
we  may  write  Eq.  (2)  as  follows t 


r Ov.yV^ 

W1'1.'  - 


0Ch(xo»R)  ^ H(*l>xo’ 


^cAl^ 


/ ln£>$ 


dx 


x 

+^K^xo*R^K^x2,xo*R^  J W dx* 


or 


<Ws 


^ OtH(xQ,R)  <Th(x1,xo)  + 0<k(xo,R)(5k(x2,xo,R), 


where,  according  to  Eq.  (A-2), 


x1  ■ J xln(^  )<§dx  j J'  ln(^)4^dx, 

o In 

x2  - J'  X $ dx  j J $ dx. 


(5) 


(6) 


(7) 


(8) 


! . 


and 


rlm  / in(j2) 


(9) 


The  three  quantities,  x^,  Xg,  and  , according  to  their  definitions,  are 

functions  of  the  residual  range  of  n -mesons  at  x ,R.  Figures  10  and  11 

° _2 

show  these  quantities  plotted  versus  R(100  <R<  2*,  000  g cm  ) for  x = 

-2  0 
1,000  g cm  . The  calculations  were  carried  out  numerically  whereas  the 

function  ^ (x,x  ,R)  was  computed  with  the  aid  of  Eqs.  (H-33)  and  (11-18) 
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which  corresponds  to  i*0°  of  geographic  latitude.  Note  that  the  factor 

_o 

is  practically  one  for  all  R > 5 00  g cm  j note  further  that  the 
atmospheric  depths  x^  end  Xg  are  slowly  deereasing  functions  of  R.  W© 
shall  discuss  the  physical  significance  of  these  results  in  Section  II-D. 

In  order  to  give  Eq.  (6)  a simpler  physical  interpretation,  it  io 
useful  to  approximate  the  function  £ K by  the  following  expressions t 


r*  x +R+b  /x  +R+b-x«\  f ° 

6k(x2,x  .R >»-|~  in  -Srl-  / CT T(x.  )dx«  • 


xo"*2 


one  may  readily  verify  that  the  above  formula  is  a good  approximation  to 
Eq.  (It* ) If  T(x')  do*s  not  vary  too  strongly  with  x'  in  the  region  be- 
tween xq  and  Xg.  In  this  oase  Eq.  (6)  simplifies  tot 


^Ti/iv  ’ aH^H(xl»xo)  + V^AV^o5  * 


where 


*h'<Vvr)V 

dK  -0^(xq,R)  (xQ+R+b)  In 4 


x +R+b-x„ 
o c 


( ^T)av  (x2,x0) 


xo_x2 


J T(x* 


Eq.  (10),  although  not  so  accurate  as  Eq.  (6),  gives  us  an  insight  into 

the  problem  in  question t the  tempera ture  effect  on  the  differential 

intensity,  i^,  at  sea  level  can  be  described  by  means  of  a two-term 

regression  formula.  The  variations  of  iy  are  to  be  correlated  with  the 

variations  of  the  height  of  and  the  variations  of  the  temperature 

averaged  between  x0  and  x . Figures  12  and  13  ahow  the  coefficients  au 
d.  o n 

and  plotted  versus  the  residual  ranges  R. 

r 


I 


* 
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We  conclude  this  section  with  a few  remarks  concerning  the  tem- 
perature effect  for  the  integral  intensity  near  sea  level,  I?.  Evi- 
dently, one  obtains  an  expression  applicable  to  this  effect  by  inte- 
grating Eq,  (6)  or  Eq,  (10)  over  all  residual  ranges  R above  certain 
cut-off  value,  Rq,  determined  by  the  experimental  arrangement,  viz t 


^tVRo,Xo)  " J a^^HdR  ♦ J 1) . 

R 


(1?) 


Since  the  functions  6*  H and  (d*T)Av  vary  very  slowly  with  respect  to  R, 
Eq,  (12)  may  be  evaluated  by  means  of  the  mean  value  theorem.  Thus,  ac- 
cording to  Eqs.  (A-l* ) and  (A-2)t 


£TV*V  " VH(£L 


,xo)  ♦ Ar(( 


T) 


) 


Av'  2 


(12*) 


where 


00 


R 

o 


^ - x1(R1)  , x2  - x2(R2)  | 


(13) 


Figures  12  and  13  show  the  coefficients  A^  and  A^  plotted  against  the 

cut-off  range,  Rq.  These  curves  are  convenient  for  comparative  purposes 

with  experimental  observations.  We  shall  return  to  them  in  Section  II-D  . 

The  numerical  values  of  5L  and  x„  for  the  case  of  I (R  ) may  be 

X <;  vo 

estimated  directly  from  Figure  10  where  one  has  to  take  those  values  of 
x^  and  x2  which  correspond  to  the  average  residual  ranges  R^  and  RgP 


I 


s 


respectively.  On*  finds  that  5L  trill  lie  somewhere  between  115  and  110  g cm  , 
— -2 

and  Xg  between  190  and  160  g cm  # depending  on  the  experimental  value  of  Rq. 

C-  Pressure  Effect  on  n-Keson  Intensity. 

To  complete  the  discussion  of  the  atmospheric  effects  on  the  cosmic-ray 
intensity  it  remains  to  evaluate  the  effect  of  the  atmospheric  pressure. 

This  evaluation  may  be  carried  out  by  a method  analogous  to  that  discussed 
in  the  preceding  section.  If  the  temperature  overhead  is  kept  constant, 
the  partial  variation  of  the  differential  intensity,  i^(R,x^),  due  to  the 
changes  of  the  ground  pressure,  xq,  is  given  by* 


r*  <^i  (R,x  ) r 

<WR’*o>  * ~^x <5*0  • 


/ 


Assuming  that  i^  may  be  represented  by  Eq.  (H-8)  one  findB  by  differentia- 
tion* 

f x /L  f°  } (lit* ) 

| ° + I ~-JinD(R')+lnw(x,xo,R^<f  dx>/ xq  . 

o 

The  first  term  in  Eq.  (lit' ) is  negligible  at  sea  level.  The  two  last  terms 
may  be  evaluated  by  means  of  the  mean  value  theorem  discussed  in  the  Appen- 
dix. By  substituting  5qs.  (H-35)  (H-lS)  for  2 (?»’ ) and  v(x-xo,R), 
respectively,  one  gets* 

fpViy  * aP(fx0  » (1^) 

where  the  "pressure*  coefficient,  *p,  is  given  by* 

•5  to  <l)0dp(xo»R) 

aP  " “ 52fieio*R^2  + ^ T&IT H(X1»X0) 

(16) 

* aH<^'R)^TSl)  * [°VVE>  «wr)J  • 

o 

The  evaluation  of  Eq.  (16)  with  the  help  of  Eqs.  (11-19)  to  (11-22)  and 
Fig.  6 shows  that  the  first  term  above  is  prominent.  We  have  computed 


1 


rl 


\?**«** 
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the  pressure  coefficient  ap  for  100  <"  R U,000  g cm  by  using  for  the 
temperature,  T(x),  the  annual  average  at  lt0°  n geographic  latitude.  The 
result  is  presented  in  Fig.  lit,  where  ap  is  plotted  versus  R» 

For  the  convenience  of  the  reader  we  have  also  computed  the  pressure 
coefficient  of  the  integral  intensity,  i.e«,  the  quantityt 

co 

AP  - j J VvdR  * (17) 


The  behavior  of  Ap,  as  a function  of  Rq,  is  shown  in  Fig.  lU. 

Before  concluding  this  section,  we  should  like  to  point  out  that 
the  coefficient  ap  is  directly  related  to  the  measurements  on  the  altitude 
dependence  of  the  differential  intensity,  ivt  according  to  its  definition 
ap  simply  represents  the  slope  of  the  intensity-depth  curve.  Therefore, 
referring  to  Eq.  (16),  the  observed  values  of  ap  may  be  alternately  used 
to  check  the  correctness  of  the  logarithmic  derivative  of  G(R')» 


D. 


Summary.  Comparison  with  Experiments. 

(a)  Atmospheric  effect  on  the  differential  intensity. 


According  to  the  results  obtained  in  Sections  B and  c the  total 
variation  of  the  differential  intensity  near  sea  level  may  be  represented 
by  a three-term  regression  formula,  vizt 


faAr  " aH^H<xi»x0)  + aK(^T)Av(Vxo)  * “p^xo  * 


(18) 


i.e.,  the  fractional  changes  in  i are  to  be  correlated  with  the  fol- 


lowing variables  of  the  atmosphere r 

1)  the  departure  from  the  mean  height,  S H,  of  the  pressure 

^(R)# 

2)  the  departures  from  the  mean  temperature  overhead  averaged 
between  the  pressure  level  Xg(R)  and  the  sea-level  pressure. 


xQ,viz: 


((f-M  x ) - t 


'AV  -2»“o'  x_-x„  J 


o 2 


T 


m 


6T(x')dx'  ; 


•cm"2  (air  equivalent 
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The  numerical  values  of  the  two  characteristic  pressure  levels,  x^  and  x2, 
are  to  be  taken  from  Fig.  10.  For  the  numerical  values  of  the  three  co- 
efficients, a^,  a^,  and  ap,  one  is  referred  to  the  Figures  12,  13  and  lU. 

(b)  Atmospheric  effect  on  the  integral  intensity. 

Similarly  to  case  (a),  the  total  variation  of  the  integral  intensity 
near  sea  level  is  given  by: 

<f  VTv  = AH<f*l(xL,xo)  + Ak(  g T)Av(x2,xo)  + Ap  <f  xo.  (19) 

Here  the  coefficients  Ap,  A^,  and  Ap  are  constant  quantities  for  a given 
experimental  arrangement.  Their  numerical  values  will  be,  however,  dif- 
ferent from  case  to  case,  depending  on  the  amount  of  shielding  material 

_9 

above  the  detector.  Taking  U00  g cm  air  equivalent  of  shielding  material 
as  a typical  case,  we  find  from  Figures  12,  13  and  li*: 

Ap  * -3.15£  per  km,  ") 

Ak  —0.0592  per  °C,  ^ (20) 

Ap=- 1.79!?  per  cm  Hg.  J 

In  order  to  determine  the  ab^ve  coefficients  experimentally,  one  has 

to  correlate  the  observed  changes  in  Iv  with  the  three  atmospheric 

variables  (ThCsc^x^),  ( <^T)Av(xp,xo)  and  cTxo.  The  characteristic 

pressure  levels,  x-(  and  x2,  have  now  practically  fixed  values:  x^-^115  g cm 

and  x„  ^190  g era  u,  if  the  amount  of  shielding  above  the  detector  is  of 
d -2 

the  order  of  a few  hundred  g cm  . It  is  worth  mentioning  that  the  pres- 
sure level  x2  nearly  coincides  with  that  of  the  tropopause  at  moderate 
latitudes.  This  implies  that  the  vertical  average  of  the  temperature 
changes  above  the  observer  ( <^*T)Av,  extends  only  over  the  region  of  the 
troposphere,  and  does  not  include  the  temperature  of  the  stratosphere. 

(c)  Comparison  with  experiments . 

As  we  have  mentioned  in  Sec.  A,  it  has  been  customary  to  correlate 
the  observed  variations  of  the  p.-meson  intensity  with  only  two  atmos- 
pheric variables,  (f*H  and  x . 


Lack  of  knowledge  of  the  third  variable, 


s 


£ T).  . 


( O T)^  , makes  an  exact  comparison  of  our  results  with  those  found 
experimentally  impossible.  However,  the  following  semi -quantitative 
considerations  appear  to  indicate  substantial  agreement  between  obser- 
vations and  theory. 

_2 

In  the  first  place,  our  value  of  115  g cm  for  the  "mean  pressure 

level*  of  n -meson  production  is  in  agreement  with  the  observations.  For 

example,  we  have  already  mentioned  that  A.  Duperier  obtained  the  best 

correlation  between  changes  in  cosmic -ray  intensity  and  changes  in  the 

height  of  a given  pressure  level  by  choosing  a value  of  this  pressure 

-2 

level  in  the  proximity  of  100  g cm  • 

In  the  second  place,  if  one  writes  Eq.  (19)  in  the  formt 

s V1, ' [v*K  sr^J  H * ap5xo  • 

one  sees  that  the  so-called  "decay*  coefficient  of  Eq.  (l),  Afo,  is 
related  to  our  coefficients  and  by  the  equation: 

(£t>Av 

Ah“  AH  * h ~$H~  * ^ 

The  ratio  (<5" 1)^/0  H will  depend,  of  course,  on  the  geographic 
location,  the  season  during  which  the  experiment  was  performed,  etc. 
Meteorological  observations  indicate  that,  for  moderate  latitudes,  Hie 
above  ratio  will  have  on  the  average  the  following  approximate  values: 
for  diurnal  departures  of  (T)^v  and  H from  their  means 

[( S T>Ay  4-h]  50  c°c/tao, 

for  seasonal  departures  of  (T).  and  H from  their  means 

r-  n AV  _ 


<5'hJ 


seasonal 


20  (°C/km). 


It  follows  that  the  decay  coefficient,  A^will  be  different  depending  on 
whether  it  is  inferred  from  the  diurnal  on  the  seasonal  changes  of  the  in- 
tensity. According  to  Eqs.  (20)  and  (22),  one  would  expect  for  A^  in  these 


two  cases: 


Ofi  ] diurnal  * 6>1  * Per  k"- 

fi]  saaaonal*-  h'3  * P*r  lan* 


U.3  % per  km. 


-hh~ 


Indeea,  the  diurnal  and  seasonal  decay  coefficients  estimated  above  are  in 
substantial  agreement  with  those  observed  experimentally.  For  instance* 

Dolbear  and  Elliot  (DDWSl)  deduced  for  A'  the  values  -5.7  % km  ^ and  -3.6  % 

-1  M 
km  from  the  daily  and  monthly  correlations,  respectively.* 

Regarding  the  pressure  effect  on  the  p-meson  intensity,  we  shall 
limit  ourselves  to  the  remark  that  the  values  of  the  pressure  coefficients 
Bp  and  V*  given  in  Fig.  ll*,  seem  to  be  in  essential  agreement  with  those  ob- 
served experimentally.  This  can  be  demonstrated  by  a comparison  of  the 

differential  range  spectra  for  two  different  depths  near  sea  level,  as, 

_2 

e.g„,  the  spectra  measured  by  Conversi  in  Chicago  (1,010  g cm  ) and  at 

_2 

China  Lake  (937  g cm  ).  One  finds  that  by  correcting  the  sea-level  spec- 
trum given  in  Fig.  2 by  means  of  ap  for  the  two  corresponding  levels  one 
obtains  two  curves  practically  identical  with  those  given  by  Conversi 
(CH50).  The  agreement  may  be  considered  as  a partial  check  for  the  cor- 
rectness of  the  production  spectrum  of  p -mesons. 

In  conclusion,  we  would  like  to  discuss  briefly  the  problem  of  the 
so-called  "positive  temperature  effect".  As  we  have  mentioned  in  Sec.  A 
recent  observations  by  Dupsrier  have  shown  that  the  two -term  regression 
formula  given  by  Eq.  (1)  was  inadequate  to  account  fully  for  the  observed 
variation  of  the  cosmic -ray  intensity.  The  analysis  of  the  correlation 
coefficients  indicated  that  there  must  be  an  additional  atmospheric 
variable  which,  together  with  H and  playB  a controlling  role  in  the 
intensity  variations.  Duperier  assumed  this  variable  to  be  the  temperature 
of  the  pressure  layer  between  100  and  200  mb  (lower  stratosphere).  Hence, 
he  replaced  Eq.  (1)  by  the  following: 

&I/I  - AjjSh  + ♦ Ap  5xo  , (23) 

where  is  the  deviation  from  the  mean  of  the  temperature  of  the 
100-200  mb  pressure  layer,  and  the  other  symbols  have  the  same  meaning 
as  before.  The  temperature  coefficient.  A™,  deduced  from  the  observational 


* In  the  experimental  arrangement  used  by  Dolbear  and  Elliot  the  threefold 
counter  coincidence  set  was  inclined  at  to  the  vertical.  Therefore, 
the  results  of  these  authors  are  to  be  corrected  before  being  compared  with 
the  results  corresponding  to  the  vertical  intensity.  It  can  be  shown,  how- 
ever, that  this  correction  is  not  crucial,  and  can  be  neglected  in  cur  semi- 
quantitative  discussion. 
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data,  turned  out  to  be  positive  (A^  ■ 0*12  % per  °C).  Duperier  attempted 
to  interpret  this  positive  temperature  effect  as  due  to  ’the  competing 
processes  of  nuclear  capture  and  decay  of  n-mesons.  However,  a quantitative 
estimate  of  the  effect  of  these  processes,  based  on  more  recent  data  for 
the  mean  life  and  the  oross  section  for  nuclear  capture  of  tt -me  sons,  »hows 
that  the  observed  value  of  is  much  too  high  to  be  ascribed  exclusively 
to  the  finite  life  span  of  w-mesons. 

It  is  interesting  to  compare  Duperier’ a regression  formula,  Eq.  (23), 
with  that  given  by  Bq.  (19).  Since  the  coefficients  and  Ap  determined 
by  Duperier  are  roughly  in  numerical  agreement  with  those  calculated  in 
this  paper,  the  ratio  ( $ T)^/  6“ T should  be  equal  to  A^/A^.  The  value 
of  our  coefficient,  A^,  is  about  -0.06  % per  °C  while  the  experimental 
value  of  Aj  is  0.12  % per  °G  from  which  it  follows  that  ( <5  T)  ,y/  -2. 

The  negative  sign  of  the  ratio  ( S T)^/  <5*T  is  not  in  contradiction  with 
the  general  behavior  of  the  free  atmosphere.  For  it  is  well  known  that 
the  warming  up  of  the  troposphere  [j  T)^v  positive]  is,  as  a rule,  ac- 
companied by  a cooling  down  of  the  lower  stratosphere  ( & T negative)  and 
vice  versa.  The  magnitude  of  the  above  ratio  is  probably  too  large.  How- 
ever, a knowledge  of  the  actual  values  of  ( S T)^v  for  the  period  and  loca- 
tion of  Duperier’s  experiment  is  needed  to  check  our  results  quantitatively. 

In  view  of  the  above  discussion  we  conclude  that  the  additional  term, 
A^(  <$*  T )^v,  in  the  regression  formula  may  well  remove  the  apparent  seasonal 
variability  of  Elliot’s  decay  coefficient  as  well  as  the  anomalous  value 
of  Duperier’s  coefficient  for  the  positive  temperature  effect.  An  ex- 
perimental verification  of  this  conclusion  would  be  desirable.  Unfortu- 
nately, because  of  the  strong  correlation  between  S' H and  ( S'  T)^.r,  it 
will  be  very  difficult  to  separate  experimentally  the  effects  caused  by 
S’ H from  those  caused  by  ( S'  T)^y. 

17.  LATITUDE  DEPENDENCE  OF  n-MESON  SPECTRUM  AT  PRODUCTION 

A.  Introductory  Remarks. 

Ever  since  the  earlisat  years  of  cosmic -ray  research  it  has  been 
known  that  the  local  charged  cosmic  radiation  changes  in  a definite  pattern 
when  one  passes  from  south  to  north  on  our  globe  (for  a review,  see,  for 
example,  KH7*>2).  This  empirical  evidence,  although  not  understood  in  all 
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detail,  ha3  had  a great  heuristic  value  in  guiding  physicists  toward  a 
more  fundamental  phenomenon » the  geomagnetic  effect  on  the  primary  cosmic 
radiation.  Today  we  are  in  a position  to  link  the  predicted  effects  of 
the  earth' s magnetic  field  with  those  observed  experimentally  and  thus  to 
strive  for  a consistent  picture  of  all  phenomena  related  to  the  latitude 
variation  of  cosmic  rays.  Among  several  striking  features  concerning  the 
behavior  of  a charged  particle  in  the  earth's  magnetic  dipole,  the  fol- 
lowing will  be  of  particular  interest  to  uei  a charged  particle  of  energy 
lower  than  a certain  limit  is  prevented  from  reaching  the  top  of  our  at- 
mosphere and  thus  is  excluded  in  the  observed  energy  spectrum  of  the  pri- 
mary radiation.  The  value  of  this  cut-off  energy  depends  on  the  charge 
and  mass  of  the  particle,  and  changes  with  geomagnetic  latitude  as  well 
as  with  the  angle  of  incidence.  Columns  2 and  3 of  Table  II  show  these 
values  for  the  case  of  vertical  protons  and  of -particles  (predominant 
constituents  of  primary  radiation).  The  observations  at  high  altitudes 
carried  out  by  means  of  rockets  and  balloons  not  only  support  t-he  theore- 
tical predictions  on  the  cut-off  effect  but  also  provide  us  with  the  ab- 
solute numbers  of  primary  particles  which  reach  the  upper  atmosphere  at 
various  geomagnetic  latitudes.  The  results  obtained  by  Winckler  et  al. 
(WJRf>0)  are  compiled  in  column  U of  Table  II. 

Since  we  are  concerned  with  the  analysis  of  the  meson  component  of 
the  cosmic  radiation,  the  problem  which  is  of  importance  to  us  is  the 
following t how  does  the  latitude  variation  of  the  primary  flux  affect,  the 
production  of  n -me sons  and  thus  the  production  spectrum  of  pi-mesons? 
Qualitatively,  it  is  evident  that  the  observed  decrease  in  the  primary 
flux  should  lead  to  a reduction  of  the  number  of  n-mesons  produced.  (The 
above  statement  would  be  untrue  only  if  primary  particles  of  energies  be- 
low the  cut-off  limits  at  the  geomagnetic  equator  were  completely  inefficient 
in  meson  production.)  However,  in  order  to  estimate  quantitatively  the 
magnitude  of  this  reduction  on  the  basis  of  the  observed  decrease  of  the 
primary  radiation,  one  would  have  to  know  the  cross  sections  for  the  pro- 
duction of  n-meson3  as  functions  of  the  energy  of  primary  particles.  Since, 
for  the  energy  region  with  which  we  arc  concernsd  here,  the  cross  sections 
are  not  known,  the  only  way  to  obtain  some  information  c the  latitude  de- 
pendence of  the  p -meson  spectrum  at  production  is  to  study  the  latitude 
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TABLE  II 

Geomagnetic  data  pertinent  to  the  primary  radi.  tion  and  the  it-  and 
p -meson  components.  The  symbols  below  have  the  following  meanings 
y)  is  the  geomagnetic  latitude;  and  are  the  kinetic  cut-off 
energies  of  vertical  primary  protons  and  a-particles,  respectively; 
j is  the  total  vertical  intensity  of  primary  radiation  (WJR50) ; 
r^  - [iv(100,300)j^  /[iv(lCO,300)j  5qo,  where  1^(100,300) 

is  the  differential  vertical  intensity  of  n -mesons  with  the  residual 

-2  -2 
range  of  100  g cm  at  the  atmospheric  depth  of  300  g cm  , for  a 

given  geomagnetic  latitude,  vj>  (CM50);  a is  the  parameter  appearing 

in  the  empirical  formula  for  the  production  spectrum  of  p -mesons : 

G(R')  ■ 7.31  x 10^  (a+R> )”3.58^  TP  is  the  total  number  of  n-mesons 

of  energies  greater  than  260  Mev  produced  throughout  the  atmosphere 

in  the  downward  cone. 
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dependence  of  measurable  quantities  which  are  related  in  a known  fashion 
to  the  p. -me  son  spectrum  at  production.  One  such  quantity  is  the  differential 
vertical  intensity  of  |x -mesons  discussed  already  in  the  preceding  sections; 
its  relation  to  the  p-meson  spectrum  at  production  is  expressed  by  Eqs.  (II-8), 
(II -18)  and  (11-33).  It  follows  from  Eq.  (II-8)  that,  from  the  known  de- 
pendence of  iv(R,s)  and  w(x,s,R)  on  geomagnetic  (or  geographic)  latitude, 
one  should  be  capable  of  drawing  some  conclusions  about  the  latitude  de- 
pendence of  the  p-meson  spectrum  at  production. 

B«  Experimental  Evidence  of  Latitude  Dependence  of  p -Meson  Intensity. 

In  contrast  to  the  abundance  of  the  experimental  material  concerning 
the  latitude  dependence  of  the  total  penetrating  cosmic  radiation  the  in- 
formation on  the  latitude  dependence  of  the  differential  p-meson  intensity 
itself  is  very  incomplete.  It  is  true  that  near  sea  level  one  can  identify 
the  total  penetrating  cosmic  radiation  with  the  p-meson  component  and  thus 
regard  it  as  usable  for  our  analysis.  However,  the  latitude  variation  is 
so  greatly  reduced  near  sea  level  as  compared  with  higher  elevations  that 
it  does  not  represent  any  promising  basis  for  quantitative  estimates  cf 
the  latitude  dependence  of  the  meson  spectrum  at  production.  As  far  as 
the  high-altitude  measurements  are  concerned,  we  know  of  only  one  experi- 
ment, namely,  that  carried  out  by  M.  Conversi  (CM^O)  aboard  an  aircraft 
in  several  series  of  flights  at  various  latitudes.  Conversi' s measure- 
ments were  made  at  the  constant  altitude  of  30*000  feet  with  the  apparatus 
and  technique  discussed  in  Part  II,  Sec.  B.  Thus  his  data  can  be  inter- 
preted as  a direct  measure  for  the  latitude  dependence  of  i (R,s)  for 
-2  -2 

R « 100  g cm  and  s ■ 307  g cm  . Unfortunately,  the  statistical  accuracy 
achieved  in  this  series  of  measurements  is  poor  (at  some  latitudes  the 
errors  exceed  20  percent)  so  that  not  much  weight  can  be  given  to  the 
shape  of  the  latitude  curve  chosen  by  Conversi  as  the  best  fit  to  the 
small  numbers  of  points  obtained.  Nevertheless,  it  is  well  to  remember 
that  this  curve  still  represents  the  best  information  on  the  latitude  de- 
pendence of  p-meson  intensity  at  high  altitudes  available  at  the  present 
time. 

Tt  will  be  useful  for  the  purpose  of  our  discussion  to  express  the 
•i-wneson  intensities  measured  at  various  latitudes  in  multiples  of  the  in- 
tensity at  50°  V geomagnetic  latitude  rather  than  at  the  geomagnetic  equa- 
tor. Referring  to  Conversi' s curve,  we  quote  these  ratios  in  column  5 of 
Table  II. 
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C,  Atmospheric  Latitude  Effect. 

It  follows  from  the  discussion  sketched  in  Sec.  A that  there  are 
actually  two  causes  for  the  observed  variation  of  the  p-raeson  intensity 
with  latitude.  While  the  first  cause  originates  in  the  geomagnetic  cut- 
off of  the  primary  spectrum  and  affects  the  p -meson  spectrum  at  production, 
the  second  cause  originates  in  the  latitude  dependence  of  the  vertical  dis- 
tribution of  the  atmospheric  temperature  and  affects  the  survival  probability 
of  p-nesens.  Although,  as  we  shall  see  later,  the  latter  effect  proves  to 
be  in  general  of  minor  importance,  in  some  special  cases  it  cannot  be  dis- 
regarded completely.  Thus,  before  attempting  any  quantitative  deductions 
concerning  the  georaagnetio  effect  on  the  production  spectrum,  we  must  assure 
ourselves  as  to  the  magnitude  of  the  second  effect. 

Let  us  first  discuss  in  general  terms  how  the  existing  temperature 
distribution  in  our  atmosphere  affects  the  intensities  of  p -mesons.  If 
we  limit  ourselves  to  the  annual  mean  conditions  on  our  globe  we  can  base 
our  considerations  on  Fig.  lx.  This  figure  implies  that  if  we  move  south- 
wards along  an  isobaric  surface  below  the  tropopause  we  encounter  in  the 
northern  hemisphere  a positive  temperature  gradient.  On  the  other  hand,  an 
excursion  above  the  tropopause  in  the  same  direction  is  confronted,  as  a 
rule,  by  a temperature  gradient  of  the  opposite  sign.  Since,  according 
to  Eq.  (11-18),  an  increase  in  values  of  the  integrals  ^ 3 T(x* )dx'/x'  and 
£ T(x*  )dx'  reduces  the  chance  of  survival  of  a p-meson,  it  follows  that 
the  existing  horizontal  temperature  gradients  can  reduce  or  enhance  the 
p-meaon  intensities  depending  on  whether  we  make  the  measurements  (moving 
southward)  below  or  above  a certain  atmospheric  depth.  The  layer  where 
this  atmospheric  latitude  effect  changes  its  sign  is  to  be  expected  some- 
where between  2$0-  and  300-mb  levels.  As  to  the  magnitude  of  the  effect, 
we  infer  further  from  Fig,  lx  that  its  absolute  value  is  greatest  for  the 
intensity  measurements  near  sea  level. 

All  the  above  findings  are  confirmed  ly  explicit  numerical  calculations. 
For  observations  near  sea  level  we  can  apply  the  regression  formula  (III-10) 
to  estimate  the  relative  changes  of  p -me son  intensities  with  respect  to  that 
at  U0°  geographic  latitude.  The  results  of  our  calculations  are  shown  in 
the  first  6 columns  of  Table  III.  One  notes  that,  whereas  for  high-energy 
mesons  the  effect  amounts  only  to  a few  percent,  the  relative  orange  of 
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table  hi 


Atmospheric  latitude  effect  on  the  p-meson  intensity  near  sea  level 
and  at  the  depth  of  300  g cm  The  symbols  PVMf  represent 
the  relative  departures  of  the  i (R,  s) -values  at  a given  geographic 
latitude,  f , from  the  corresponding  i^(R, s) -values  at  U0°  geographic 
latitude . 


IS i (R,x  )/i  (R,x  ) I 
L v'  1 o ' v'  ’ o J||)i 


5oo 

■| 

1,000 

1,500 

2,000 

-o.Q3lil 

-O.0196 

-0.0132 

-0.0093 

-0.0155 

-0.0089 

-0.0059 

-O.OOU2 

0 

0 

0 

0 

0.0156 

0.0087 

o.oo56 

0.0038 

0.0283 

0.0170 

0.0119 

0.0088 

(100,303) 


iv(100,300)j  . 

-0.016U 

-0.0055 
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0.0057 

o.oU*5 
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slow-meson  intensities  can  b#  considerable.  They  are  certainly  not 
negligible  if  one  compares  the  differential  intensities  of  ^-mesons 
with  R ■ 100  g cm  ‘ observed  at  middle  latitudes  with  those  observed 
near  the  equator. 

Regarding  the  atmospheric  latitude  effect  at  high  elevations  we  are 

particularly  interested  in  the  quantitative  estimate  of  Toi ^/i  , 

-2  -2  u » v 

at  s ■ 300  g cm  and  R ■ 100  g cm  (Conversi's  measurements).  Since 

for  this  atmospheric  depth  the  regression  formula  (III-10)  is  not  applicable 
[see  condition  (2)  in  the  Appendix]  we  have  carried  out  rigorous  numeri- 
cal computations  of  iv(10G,  300)  for  various  geographic  latitudes  by  making 
use  oi  Eq.  (II-8)  and  Fig.  U.  The  results  expressed  in  terras  of  relative 
departures  from  the  iy(l00,  300) -values  at  U0°  geographic  latitude  are 
given  in  the  last  column  of  Table  III.  One  sees  that  the  effect  is  small, 
as  expected  from  the  above  qualitative  considerations. 

Before  closing  this  section  we  would  like  to  clarify  one  point t in 
our  numerical  calculations  we  have  tacitly  assumed  that  the  production 
spectrum  G(R' ),  is  identical  at  all  latitudes  and  is  given  by  Eq.  (11-33). 
Although  this  assumption  as  such  is  certainly  wrong,  it  can  be  shown  that 
the  quantities  quoted  in  Table  III  are  practically  unaffected  by  the  vari- 
ation of  Q(R'  ) . Consider,  for  instance,  the  sea-level  data  computed  on 
the  basis  of  Eq.  (III-10).  Among  the  quantities  appearing  on  the  right- 
hand  side  of  this  latter  equation,  only  ^ , x^,  and  Xg  are  affected  by 
the  assumed  values  of  0(R'  ).  If  one  takes  into  account  their  insensitivity 
with  respect  to  G(R')»  and  the  manner  in  which  they  enter  Eq.  (IH-10), 
one  concludes  that  the  relative  intensity  changes  caused  by  horizontal, 
temperature  gradients  at  constant  pressure  are  practically  independent  of 
the  actual  variations  of  the  production  spectrum. 


D.  Geomagnetic  Effect  on  the  p -Meson  Spectrum  at  Production. 

The  discussion  of  the  foregoing  sections  indicates  that  there  are 
several  advantages  which  are  offered  by  the  study  of  the  latitude  dependence 
of  p.-meson  intensities  at  high  altitudes.  Among  others,  the  following  are 
worth  mentioning*  (1)  the  effect  is  sufficiently  large  to  be  recognized 
as  real  even  with  considerable  statistical  and  instrumental  errors:  (2)  the 

atmospheric  latitude  effect  is  insignificant.  In  other  words,  the  latitude 
dependence  of  p. -meson  intensity  in  the  proximity  of  300  mb  originates 


almost  exclusively  from  the  geomagnetic  effect  on  the  n-meson  spectrum 

at  production.  We  may,  therefore,  use  Conversi's  latitude  curve  as  a 

direct  index  for  the  latitude  variations  of  this  quantity.  Although  it 

is  evident  that  Conversi's  curve  alone  is  not  sufficient  to  determine 

completely  the  behavior  of  G(R' ) at  all  latitudes,  it  can  be  of  great 

value  if  it  is  supplemented  by  s?ome  additional  information.  One  such 

information  is  the  fact  that  at  large  residual  ranges  the  production 

spectrum  is  nearly  independent  of  the  geomagnetic  latitud i.  This  follows 

from  the  nature  of  the  geomagnetic  effect.  A (i-meson  wit):  a residual 

-2 

range  larger  than,  say,  3*000  g cm  has  for  its  parent  a n-meson  of 

kinetic  energy  of  not  less  than  8 Bev  (see  Fig.  16).  A n-meson  of 

that  high  energy  can  be  produced  by  the  primary  particle  only  if  the 

energy  of  the  latter  exceeds  12  Bev.  Thus,  most  of  the  ^-mesons  with 

-2 

residual  ranges  above  3 >000  g cm  are  originated  by  primaries  which 

are  sufficiently  energetic  to  overcome  the  magnetic  field  of  the  earth 

at  all  latitudes.  Hence,  the  production  spectrum  of  ^-mesons  shown  in 
-2  * -2 

Fig.  8 (3,000  g cm  <f  R'  < 7.000  g cm  ) represent  a solution  not  only 
for  50°  geomagnetic  latitude  but  for  all  latitudes. 

Let  us  make  a tentative  assumption  that  the  production  spectrum, 

G(R' ),  can  be  expressed  at  all  latitudes  by  a formula  of  the  type  given 
by  Eq.  (11-32).  Then  we  see  that  among  the  three  characteristic  con- 


stants, C.  a and  y , only  a can  depend  on  the  geomagnetic  latitude  as  a 
is  the  only  constant  with  respect  to  which  G(R')  is  insensitive  at  large 
ranges.  Expressing  this  in  symbols  we  expect  G(R')  to  be  given  by: 


h 

- S- 


G(R' ) 


; 7.31  x 1CT 
[a(tp  )-*-Rj3,! 
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I cm  sec  sterad 


where  indicates  the  geomagnetic  latitude.  The  determination  of  the 
latitude  dependence  of  a(  <p  ) represents  a relatir~l^  simple  problem. 

By  computing  the  right-hand  side  of  Eq,  (II-8)  for  various  values  of 
a(  ip  ) at  constant  R and  s and  by  comparing  the  results  with  the  i^CR, s)  - 
values  pertinent  to  given  geomagnetic  latitude,  £i^(R,s)J  ^ , one  arrives 
at  the  functional  relation  between  a(  y?  ) and  [±v(R,s )J^  . Referring  to 
Conversi's  latitude  curve  we  show  in  Fig,  15  the  functional  relation  between 
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a(  ) and  r^>  where  r^>  is  the  ratio  of  [iv(100,  300)]  ^ to 
[iy(lOO,  300)J  ^qO,  corrected  for  the  atmospheric  latitude  effect 
% (see  column  6 in  Table  IH).  With  the  help  of  this  figure  and  column  5 

of  Table  II  one  then  finds  e(  ) as  a function  of  (see  column  6 in 
Table  II).  One  notes  that  a(  ij)  ) varies  from  6U6  g cm  “ at  the  geomag- 
netic equator  to  513  g cm  ^ at  60°  geomagnetic  latitude.  Among  all  these 
values  of  a,  the  value  of  520  g cm  c at  50°  geomagnetic  latitude  is  to  be 
considered  as  most  reliable.  Its  correctness  has  been  verified  in  Part  II 
where  we  had  at  our  disposal  a relatively  accurate  and  complete  set  of 
experimental  data. 

It  is  desirable  to  test  how  plausible  is  our  assumption  that  G(R' ) 
may  be  expressed  at  all  latitudes  by  Sq.  (1).  This  can  be  done,  for 
example,  by  comparing  the  iy -values  observed  at  depths  s different  from 
300  g cm  ^ and  latitudes  <f>  different  from  50°  with  corresponding  i^-values 
! computed  on  the  basis  of  Eqs.  (l)  and  (II-8).  L, del  Rossario  and  J.  Davila- 

1 Aponte  (RL52)  have  measured  the  range  spectrum  of  slow  mesons  at  sea  level 

, at  29°  N geomagnetic  latitude.  They  found  that  i (100,x  )=(U.25  ♦ 0.13)  x 

) VO  0 

f 1(j  (g-see-sterad)  . Our  calculations  for  the  same  latitude  a(29  ) - 

_2 

595  g cm  yielded,  after  the  correction  for  the  atmospheric  latitude 

-6  -1 
effect,  iy(100,xo)  ■ U.17  x 10  (g-sec-steard)  . 

V.  PRODUCTION  SPECTRUM  OF  tr -MESONS 

i 

j A.  Significance  of  n-Meson  Spectrum.. 

! The  topics  discussed  thus  far  have  emphasized  primarily  the  "prac- 

I 

; tical"  usefulness  of  the  (i-meson  spectrum  at  production,  Q(R')«  With  the 

aid  of  G(R' ) we  were  able  to  predict  the  numerical  values  of  quantities 
which  could  be  measured  directly  by  existing  experimental  techniques. 

Now  we  turn  to  considerations  which  point  in  a different  direction.  With 
the  aid  of  G(R' ) we  shall  attempt  to  derive  a quantity  which  is  not  accessible 
> to  present  cosmic-ray  experiment  but  has  a direct  physical  meaning!  this 

l 

i quantity  is  the  energy  spectrum  of  charged  Tf-mesons  at  production.  Its 

importance  stems  from  several  reasons..  To  mention  a few,  the  knowledge 
of  the  energy  spectrum  of  charged  n-mesons  throws  more  light  upon  the 
propel ties  of  high-energy  nuclear  interactions*  by  analogy  arguments  it 
leads  tc  the  production  spectrum  of  neutral  n-mesons  and,  thus,  to  the 
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source  function  for  the  electronic  component.  These  problems  are  actually 
beyond  the  scope  of  this  analysis.  We  should  like*  however,  to  touch  at 
least  upon  one  topic*  the  multiplicity  problem  of  rr-meson  production  from 
nuclear  interactions  in  the  energy  region  of  several  Bev. 

B.  Evaluation  of  Energy  Spectrum  of  rr-Mesons  from  the  Range  Spectrum 
of  p -Mesons  at  Production. 

Owing  to  the  experimental  evidence  that  the  rr-meson  decay  is  a two- 
body  process  (tp^  -» p£  ♦ V ) we  can  derive  a relation  between  the  energy 
spectrum  of  the  p -meson  at  decay,  PCu^),  and  the  energy  spectrum  of 
p-mesons  at  production,  M(U)-  We  have  quoted  this  relation  already  in 
Part  II,  Sec.  B,  Eq.  (II-3)»  Since  in  our  case  P(U^)  is  the  unknown 
quantity,  it  is  desirable  to  solve  Eq.  (II-3)  with  respect  to  P(u  ). 

IT 

According  to  Ascoli  (AG50)  this  solution  can  be  represented  in  terms  of 
the  following  infinite  aeries t 


where 


D(U) 


p(utr)  - DOJ,)  ♦ d(u2)  ♦ ...  , 

(m  m \ 

rap  V 


r ^ , 


a) 

(2) 


and  U^,  Ug,  i<  a sequence  of  p-meson  energies  uniquely  determined 

by  the  tr-meson  energy  Ufl.  The  numerical  values  of  U^,  U2,  U,  ...  can  be 

obtained  from  Fig.  16  by  the  following  procedure t 

Beginning  with  a vertical  line  corresponding  to  a given  value  of  U^, 

find  which  represents  the  highest  energy  of  a p-moson  born  from  a 

tr-meson  of  energy  U . Proceeding  from  this  point  with  straight  lines, 
rr 

alternately  horizontal  and  vertical,  inscribe  a step  curve  within  the 
smooth  curve  drawn  in  Fig.  16.  The  points  obtained  on  the  upper  branch 
of  the  smooth  curve  then  represent,  in  increasing  order,  the  required 
sequence  U^,  Ug,  ... 

In  order  to  obtain  an  explicit  expression  for  Eq.  (2)  we  recall  that, 
according  to  Eq.  (II-5),  M(U)  - S(R')/k(U)j  thus* 

MU) 1 [Q(R„  Mil)  _dG(Rl)’] 

du  k^)  L dU  dR’  J' 


(3) 
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figure  16 

Saergy  limit*  of  a u-aeeon  born  from  e'J£-me*on  of  energy 
Un>.  TJ_and  are  the  total  energies  measured  in  multiple# 
0 i aMce  and  in'  c^,  respectively. 
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where  the  residual  range,  R',  is  now  to  be  considered  as  a function  of  U. 
The  derivative  of  the  collision  loss,  k(u),  with  respect  to  U was  com- 
puted from  the  theoretical  expression  for  k(U)  QlB52,  page  2ii,  Eq.  (7)J  ♦ 
We  found  by  differentiation* 


1 dk 

k HU  ‘ 
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A closer  inspection  of  Eq.  (3)  shows  that  the  function  D(U)  decreases 
rapidly  for  all  energies,  U,  larger  than  a certain  value  so  that  the 
series  given  by  Eq.  (1)  converges  fairly  rapidly  and  is,  thus,  convenient 
for  the  numerical  evaluation  of  P(U  ).  We  have  carried  out  the  numerical 

it 

calculations  concerning  the  energy  spectrum  of  tt -mesons  at  decay,  P(U  ), 
by  making  use  of  Eq.  (IV-1),  i.e.,  we  have  derived  P(U  ) not  only  as  a 
function  of  U but  also  as  a function  of  the  geomagnetic  latitude.  The 

tf 

results  are  summarized  in  Table  IV.  It  is  evident  from  the  previous  con- 
siderations that  the  values  of  P(Ufl)  corresponding  to  £0°  geomagnetic 
latitude  are  to  be  looked  upon  as  most  reliable.  It  is,  therefore,  worth- 
while to  study  the  spectrum  at  this  latitude  in  more  detail.  Fig.  17  shows 
the  quantity  S^(E)  plotted  versus  the  kinetic  energy,  E,  measured  in  Mev. 
This  quantity  la  related  to  P(Ufl)  simply  by* 


r<y 


(5) 


The  curve  labeled  by  S^(E)  in  Fig.  17  represents  the  corresponding 

energy  spectrum  of  n-mesons  at  production.  The  dashed  parts  of  the 

curves  in  the  low-energy  region  were  obtained  by  extrapolating  Eq.  (11-33) 
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down  to  about  10  g cm  • It  is  interesting  to  note  that  the  S^CE)  displays 
a possible  maximum  at  E “ 95  Mev.  Another  fact  which  is  worth  mentioning 
is  that  the  energy  spectrum  of  n-mesons  presented  in  Fig.  17  behaves  in  a 
maimer  very  similar  to  that  observed  at  high  altitudes  in  photographic 
emulsions  (CU50). 
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C.  Remarks  on  Multiplicity  Problem. 

The  n-meson  spectrum  at  decay  derived  in  the  foregoing  section  is 
almost  identical  to  that  at  production!  this  is  at  leaat  true  for  the 
energy  region  for  which  we  have  established  the  validity  of  Bq.  (lI~8)o 
Thus,  by  making  use  of  Table  17  we  can  estimate  the  total  number  (per 
cm  -sec-sterad)  of  tr-mescns  produced  by  the  N-component  at  given  geo- 
magnetic latitude.  Of  course,  since  we  do  not  know  the  behavior  of  tne 
production  spectra  of  n-mesons  at  energies  smaller  than  a certain  value 
Eo  we  are  capable  of  computing  only  the  following  quantity: 

CO 

7Tog  - 1 f Sff(E)dE  (Eq£,  260  Mev).  (6) 


According  to  Eq.  (6)  and  the  definition  of  S^E),  T T (Eq)  represents 
the  total  number  of  charged  rr-mesons  with  energies  greater  than  Eq 
produced  in  the  downward  cone  throughout  the  atmosphere  at  given  geo- 
magnetic latitude,  <p  . The  numerical  values  of  IT  (Eq)  for.  Ey--*  26©'  Mev 
ara  given  in  the  last  column  of  Table  II.  By  comparing  this  column  with 
column  U of  Table  II  we  note  that  the  rate  with  which  TT(  Eq)  increases 
with  geomagnetic  latitude  is  much  smaller  than  that  of  the  integral  in- 
tensity, j,  of  the  primary  radiation.  This  implies  qualitatively  that, 
on  the  average,  the  high-energy  primaries  (E^  > 1 6 Bev)  produce  mesons 
with  a higher  multiplicity  than  the  primaries  of  energies  between  2 and 
16  Bev. 

One  can  get  a more  quantitative  idea  as  to  the  dependence  of  n-meson 
production  upon  the  energy  of  the  primary  particle  by  means  of  the  fol- 
lowing argument: 

Assume  that  the  contribution  of  the  secondary  N-particles  (neutrons, 
secondary  protons  .and  n-mesons)  to  the  production  of  n-mesons  with  energies 
greater  than  Eq  is  negligible  in  comparison  with  that  of  the  primary  radia- 
tion. Then  the  quantity: 


fi  (E  ,E') 
o’  p 


aTT 

a i 


(7) 


may  be  interpreted  as  the  number  of  charged  n-mesons  with  energies  greater 
than  Eq  produced  in  the  downward  cone  by  one  primary  particle  with  energy 


s 
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where  is  ths  effective  cut-off  energy  of  the  primary  flux,  j,  at 


given  geomagnetic  latitude.  Strictly  speaking,  represents  a properly 


weighted  average  of  the  cut-off  energies  of  individual  constituents  of 
primary  radiation  (protons,  Ot-particles,  etc.).  The  quantitative  data 
concerning  the  composition  of  primaries  are  as  yet  very  wanting.  However, 
we  may  assume  with  some  certainty  that  most  of  the  mesons  will  be  produced 


by  protons  so  that  the  values  of  E^  should  not  depart  very  markedly  from 

the  cut-off  energies  of  protons,  E^,  quoted  in  column  2 of  Table  II.  Thus, 

by  identifying  E'  with  E , we  can  evaluate  Eq.  (?)  for  E ■ 260  Mev  as  a 
P P o 

function  of  E’  if  we  make  use  of  columns  7,  U and  2 in  Table  II.  We  have 

estimated  the  values  of  (All  / A j)  by  a graphical  differentiation  of 


the  curve  representing  TT  as  a function  of  j.  The  results  are  shown  in 
the  first  column  of  Table  V.  Of  course,  due  to  the  large  errors  involved 
in  both  TT  and  j,  the  multiplicity  values  quoted  in  Table  V cannot  be  given 
much  weight.  (The  simplifying  assumptions  made  above  reduce  further  the 
reliability  of  these  values.)  Nevertheless,  we  believe  that  our  results 
disclose  some  interesting  features  concerning  the  energy  dependence  of  the 


multiplicity  n(EQ,E^).  In  particular,  one  sees  that  the  increase  of  n with 


increasing  E^  is  much  more  rapid  at  lower  than  at  higher  energies.  The 


flattening-off  begins  at  about  6 Bev.  It  is  interesting  to  compare  these 
results  with  those  predicted  by  Fermi’s  statistical  theory  of  meson  pro- 
duction in  nucleon-nucleon  collisions  (FE$0).  According  to  this  theory, 
the  probability  to  observe  n tt -me sons  (both  charged  and  neutral)  produced 
lucleon-nucleon  collision  is  given  by: 


f (w)  « 
n' 


A 


t«w] n 

>i  L “ J 


(8) 


where 


E 


w 


2 y'L  + — 

» 2Mc 


is  the  total  energy  carried  by  both  nucleons  before  the  collision  in  the 

center-of -mass  system  (w  is  expressed  here  in  terms  of  the  rest  energy  of 
9 

the  nucleon,  Kc")  and  A is  the  normalization  constant  given  by: 


oo 
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(9) 
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TABLE  V 

Multiplicity  of  it -mesons  produced  by  primaries  of  kinetic  energies  of  sev- 
eral Bev$  n is  the  empirical  average  multiplicity,  i.e.,  the  total  number 
of  n-mesons  with  energies  greater  than  260  Mev,  produced  throughout  the 
atmosphere  by  a primary  particle  of  energy  E^,  and  emitted  in  the  down- 
ward direction^  n^  represents  the  theoretical  multiplicity  predicted  by 
Fermi's  statistical  model  of  nucleon-nucleon  collisions. 


E» 
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Bev 

n 

nt 

2 

0.1 

0.36 

h 

0.5 

0.66 

6 

l.l 

0.95 

8 

1.3 

1.03 

10 

1.5 

1.18 

12 

1.6 

1.30 
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By  means  of  Eqs.  (8)  and  (9)  one  oan  compute  the  expeoted  average  number 
oi  charged  n -mesons  produced  in  the  foreward  oone  in  the  nucleon-nucleon 
collision,  viz.t 


<yO 


n. 


t2-j  nfn(w)  . 

" n-0  n 


(10) 


The  factor  1/3  in  Eq.  (10)  arises  from  the  fact  that  only  2/3  of  all 
tr-mesons  are  charged,  1/2  of  which  are  contained  in  the  foreward  cone 
of  production*  Using  Eq.  (10)  we  have  computed  the  theoretical  multiplicity 
n^  for  various  primary  energies  of  interest|  the  results  are  presented  in 
the  second  column  of  Table  V.  One  sees  that,  at  lower  energies,  n^  is  lar- 
ger by  about  a factor  of  as  compared  with  our  empirical  multiplicity, 
n.  At  higher  energies  (U  Bev  Ep  <C13  Bev)  the  discrepancy  is  less  pro- 
nounced. 

Realizing  the  whole  crudeness  *rith  which  both  multiplicities,  n^  and 
n,  are  derived  we  do  not  attempt  to  draw  any  definite  conclusions  with  re- 
gard to  the  obtained  discrepancy  at  lower  energies.  However,  we  would 
like  lo  remark  that  nfc  should  be  larger  than  ii  at  lower  energies  due  to 
the  fact  . that  n contains  only  ir-mesons  with  energies  greater  than  Eo=260  Mev. 
Furthermore,  it  is  likely  that  Fermi's  multiplicity  will  be  reduced  at 
lower  energies  if  one  takes  into  account  the  plural  processes  which  occur 
inside  the  air  nucleus  during  the  nucleon-nucleon  collision.  Such  a re- 
duction is  expected,  for  instance,  from  a simple  model  of  the  air  nucleus 
suggested  recently  by  U.  Haber-Schaim  (HS5l).  Finally,  the  empirical 
average  multiplicity  contains  contributions  of  n -mesons  stemming  from  the 
second,  third  and  higher  generations.  These  contributions  oan  be  considered 
as  negligible  only  at  lowest  energies  of  primary  particles*  at  primary  ener- 
gise greater  than  10  Bev  they  become  as  important  as  those  stemming  from  the 
first  generation. 


APPENDIX 


Consider  the  following  definite  integral: 

b 

F(a,b)  - J g(z)f(z)dz  , 


(A-l) 


where  g(z)  and  f(z)  are  bounded  functions  in  the  closed  interval  (a,b). 
According  to  the  mean-value  theorem  of  integral  calculus  one  may  write 
for  F(a,b)  also:  , 

D 

F(a,b)  = g(  J )J  f (z)dz  , (A-l') 

a 

where  ^ is  a specific  value  of  z lying  somewhere  between  a and  b.  In 
general  the  value  of  ^ is  not  better  known.  However,  in  some  special 
cases  one  can  estimate  the  value  of  J fairly  accurately.  This  is  so, 
for  example,  if  the  function  g(z)  and  f(z)  satisfy  the  following  con- 
ditions, respectively: 

(1)  g(z)  is  a slowly  varying  analytical  function  in  the 
closed  interval  (a,b): 

(2)  f(z)  displays  a single  sharp  maximum  in  the  closed  inter- 
val (a,b). 

In  this  case  it  is  convenient  to  expand  the  function  g(z)  in  the  integrand 
of  Bq.  (A-l)  into  Taylor's  series  at  the  point  , determined  by  the  fol- 
lowing equation: 


•J  w 

Sj  f(z(dz  — J zf(z)dz 


(A-2) 


At  this  point  Eq.  (A-l)  becomes: 

b 


F(a,b)  - /[g(  ) ♦ \ (z-  |')2g'»(  $')♦  ...  J f(z)dz  . (A-3) 


Now,  owing  to  condition  (1),  one  may  neglect  to  a good  approximation  the 
terms  involving  the  derivatives  higher  than  the  first.  By  doing  so  and  by 
setting  one  sees  that  Eq.  (A-3)  becomes  identical  to  Eq..  (A-l'). 
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Thfc  relative  error  which  one  makes  by  neglecting  higher  derivatiM es 
of  g(z;  is  of  the  order  of  magnitude  of 


relative  error 


1 


E"  ( S ) 

■grp 


[6S ' f] 


(A-U) 


where 


Q2 J f(z)dz  - J 

a a 


r 2 

f(z)dz  » ( z f(z)dz  . 


Owing  to  condition  (2)  the  value  of  (5*  will  not  differ  markedly  from 
that  of  £ , Thus,  the  relative  error  can  be  reduced  considerably  if 
the  two  conditions  above  are  well  satisfied. 
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